Probing multiband superconductivity by point-contact spectroscopy 
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Point-contact spectroscopy was originally developed for the determination of the electron-phonon 
spectral function in normal metals. However, in the past 20 years it has become an important 
tool in the investigation of superconductors. As a matter of fact, point contacts between a normal 
metal and a superconductor can provide information on the amplitude and symmetry of the energy 
gap that, in the superconducting state, opens up at the Fermi level. In this paper we review 
the experimental and theoretical aspects of point-contact spectroscopy in superconductors, and we 
give an experimental survey of the most recent applications of this technique to anisotropic and 
multiband superconductors. 

PACS numbers: 



I. INTRODUCTION: POINT-CONTACT 
SPECTROSCOPY 

Point-contact spectroscopy (PCS) was developed more 
than 35 years ago as an experimental tool to investi- 
gate the interaction mechanisms between electrons and 
phonons in metals. Yanson^ was the first to observe that 
small microconstrictions between two metals show non- 
linearities in the I-V characteristic (and in the second 
derivative (PV/dP) that are the hallmark of inelastic 
scattering of electrons by phonons. The point-contact 
technique was later used to study all kinds of scatter- 
ing of electrons by elementary excitation in metals, like 
magnons and so on-^^. When one of the sides of a 
point contact is a superconductor, quantum phenomena 
such as quasiparticle tunneling or Andreev reflection (see 
Sect llVi^ occur at the interface, depending on the height 
of the potential barrier between the two electrodes. As a 
result, the I — V shows - in addition to the features re- 
lated to inelastic electron scattering - much stronger non- 
linearities that give rise to particular structures in the 
first derivative dl/dV (that is, in the differential conduc- 
tance) which contain fundamental information on the ex- 
citation spectrum of the quasiparticles, i.e. on the super- 
conducting energy gap and its properties in the direct and 
reciprocal space. For this reason, and apparently in spite 
of its simplicity, point-contact spectroscopy has become 
an important, sometimes unique, tool for the investiga- 
tion of superconducting materials. In some recent cases, 
PCS has provided precious spectroscopic information on 
newly discovered superconductors when more complex, 
technologically demanding techniques such as scanning 
tunneling microscopy (STM) and angle-resolved photoe- 
mission spectroscopy (ARPES) were still hindered by the 
absence of single-crystal samples of sufficient size. There 
is a number of excellent reviews that deal with the the- 
oretical and experimental aspects of point-contact spec- 
troscopy in normal metals and superconductors2i2ii^. An 
extensive and comprehensive review was dedicated espe- 
cially to point-contact results in cuprates-. The present 
review is therefore focused on the most recent applica- 
tions of point contact spectroscopy to the study of multi- 
band superconductors. A general theoretical introduc- 



tion is provided, whose aim is to explain in a simple, 
experimental-oriented way, and with a consistent nota- 
tion, theoretical models of increasing complexity for the 
interpretation of point-contact data in superconductors. 



II. FABRICATION OF POINT CONTACTS 

A point contact is simply a contact between two met- 
als, or a metal and a superconductor, whose radius is 
smaller than the electron mean free path, and this in 
most cases means that the contact is nanometric. His- 
torically, point contacts were fabricated in a number of 
ways^. The pioneering technique exploited by Yansoni 
for PCS was based on the realization of microshorts in 
the dielectric layer of a tunnel junction between two met- 
als. Another technique widely used especially in super- 
conductors (but that allows only the creation of homo- 
contacts between two electrodes of the same material) 
is the break-junction technique in which a single sam- 
ple is broken at low temperature into two pieces that 
are then brought back in contact. More recently, point 
contacts have been made by lithographical creation of a 
small hole in a thin membrane on both sides of which a 
metal film is then deposited. But the most used tech- 
nique simply consists in bringing the two electrodes in 
contact by using a micromechanical apparatus. In the 
most common configuration, often called "needle-anvil" , 
the sample to be studied is one of the electrodes, and the 
other is a metallic tip, electrochemically or mechanically 
sharpened, which is gently pressed against the sample 
surface (figure [T^) . Typically, the tip has an ending di- 
ameter of some tens of micrometers and it is easily de- 
formed during the contact^. This means that, except in 
very special cases^, parallel contacts are very likely to 
form between sample and tip^. In general this is not 
detrimental to spectroscopy, unless the sample is highly 
inhomogeneous on a length scale comparable with the tip 
end^. The needle-anvil technique has several advantages: 
i) it is non-destructive and several measurements can be 
carried out in the same samples; ii) the resistance of the 
contact can be controlled to some extent by fine tuning of 
the pressure applied by the tip. Its main drawbacks are 
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the poor thermal and mechanical stability of the junc- 
tion and the fact that, if the sample is very small (tens of 
micrometers, as it can happen with single crystals), the 
whole procedure becomes extremely difficult. For these 
reasons, since 2001 we adopted the so-called "soft" point- 
contact technique, in which the contact is made between 
the clean sample surface and a small drop (about 50 /Ltm 
in diameter) of Ag paste or a small In flake. The Ag 
or In counterelectrode is connected to current and volt- 
age leads through a thin Au wire (10 - 25 /xm in diam- 
eter) stretched over the sample, as depicted in Fig[Hb). 
Despite the large "footprint" of the counterelectrode (in 
particular in the case of Ag paste) if compared to the 
electronic mean free path, these contacts very often pro- 
vide spectroscopic information. This clearly means that, 
on a microscopic scale, the real electrical contact occurs 
only here and there through parallel nanometric chan- 
nels connecting the sample surface with the In flake or 
with individual grains in the Ag paste, whose size is 2- 
10 fj,m. With respect to the needle-anvil technique, the 
"soft" one does not involve any pressure applied to the 
sample and this can be sometimes very useful, as we will 
show in Sect. IVII CI The resistance of the as-made con- 
tacts is usually already in the suitable range for Andreev 
reflection to occur. If needed, it can be tuned by ap- 
plying short (« 50 ms) voltage or current pulses until 
a spectroscopic contact is achieved. This effect (some- 
times called "fritting"—) is well known in standard elec- 
trotechnics. The pulses have the effect of destroying some 
of the existing microj unctions and/or creating new ones 
by piercing a small oxide layer on the surface of either 
electrode. The contacts are mechanically and thermally 
very stable so that, for example, PCS measurements can 
be performed even in a cryocooler. Moreover, they can 
be made also on the thin side of small single crystals 
allowing directional point-contact spectroscopy even in 
samples too small for the needle-anvil technique. Often 
(but not always) the conductance curves of "soft" point 
contacts are more broadened than those obtained by the 
needle-anvil technique. As we will show, this is probably 
related to inelastic scattering near the interface, possi- 
bly by an oxide layer on the surface of Ag grains or of 
the sample. As a matter of fact, the same holds for con- 
tacts made with the Au wire alone, or even with a tip, 
whenever the pressure applied by the tip on the sample 
is small. 



III. POINT-CONTACT SPECTROSCOPY (PCS) 
IN THE NORMAL STATE 
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FIG. 1; (a) Experimental arrangement for point-contact mea- 
surements with the typical needle-anvil technique. A metallic 
tip (Au, Pt-Ir, Pt, Ag) is gently pressed against the surface 
of the sample, (b) The "soft" point-contact technique where 
a tiny spot of Ag paste (or a tiny flake of In) replaces the tip. 



vant quantity is the Knudsen ratio K — £/a, where i is 
the electron mean free path {!/£ = l/£e + where ie,i 
are the elastic and inelastic mean free paths) and a is the 
contact radius. From now on it will be assumed that the 
shape of the contact is a circular orifice with radius a in 
an otherwise completely reflecting barrier. Unless oth- 
erwise specifled, we will specially refer to homocontacts 
(i.e. contacts between two electrodes made of the same 
metal). Depending on the value of the Knudsen ratio, 
different regimes of conduction are possible, as described 
in the following. 



1. Ballistic regime 

In the ballistic regime the electron mean free path £ is 
much larger than the contact radius a {£ ^ aor K ^ 1). 
The applied voltage V accelerates electrons within the 
distance of a mean free path. The electrons will then 
flow through the contact ballistically (with no scatter- 
ing) gaining a kinetic energy equal to eV (see Fig. [2] 
(a)). In this way, the energy of the injected electrons is 
perfectly known and corresponds to the voltage applied 
to the junction. The resistance of the contact in such a 
situation was calculated by Sharvinl*^ and is equal to 
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where p is the resistivity of the material under study. 
Since in metals p a: £~^, Rg is independent of the elec- 
tron mean free path, and depends only on the contact 
geometry. As a matter of fact, it can also be written as 
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The uniqueness of point-contact spectroscopy in the 
normal state is due to its ability to provide spectroscopic, 
energy-resolved information on the inelastic scattering of 
quasiparticles with elementary excitations like phonons, 
magnons and so on by using a very simple and cheap 
experimental setup. To do so, however, some important 
experimental requirements must be fulfilled. The rele- 



being kp the Fermi momentum^ ^. In the k space, the 
(supposed spherical) Fermi surface (FS) expands for for- 
ward electrons by a quantity eV (see Figl^fa)). Inelastic 
scattering events taking place in the bottom electrode 
give rise to a measurable (negative) corrections to the 
current only if they cause the backflow of carriers through 
the orifice. The backscattered electron must jump back 
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FIG. 2: Schematic diagram of tlie contributions to the cur- 
rent in a point contact obtained by solving the Boltzmann 
equation, (a): zeroth-order Sharvin current (no scattering), 
(b): Diffusive regime (only elastic scattering in the contact 
area), (c): Thermal regime with inelastic scattering in the 
contact region. The right sides of the figures show electron 
distribution functions at the center of the contact for the three 
main regimes, (a): ballistic regime. The FS is formed by two 
half-spheres with different radius, i.e. defined by the surfaces 
at energy E and E + eV. (b): diffusive regime. The elastic 
scattering redistributes the electrons over the sphere but only 
in an energy shell with a width given by eV. (c): thermal 
regime. The inelastic scattering reduces the shift in energy 
space as it is usual for normal transport in a conductor. 



onto the shrunk FS and this is possible only if it can lose 
an energy eV in the scattering process. This explains 
why in the ballistic regime the applied voltage sets the 
energy scale of the spectroscopic investigation. 

The first-order correction to the current due to the 
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neffN{0) / dE / dE'S{E-E') (3) 



where ^eff = 8a'^/3 is the effective volume in which the 
inelastic scattering of electrons contributing to 51 occurs, 
A^(0) is the density of states at the Fermi level and 
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(4) 

is the spectral function for the relevant interaction, which 
results from an integration over all the initial and final 
electron states of the scattering matrix elements Iffj^j^'l 
weighted by an efficiency function Kiyi^) which ac- 
counts for the direction of the incoming and the inelasti- 
cally scattered electron. It can be shown^'^ that 
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A direct determination of the spectral function by means 
of PCS I — V measurements is thus possible. If 
the elementary excitations are phonons, S{eV) is the 
so called "point-contact electron-phonon spectral func- 
tion" a'^p(jF{eV) which differs only slightly (due to the 
efficiency function A'(k, k')) from the thermodynamic 
Eliashberg function a^FleV). In this case, using the 
formulas of the free electron model, one obtains: 
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It is worth mentioning that, according to eqlH one ex- 
pects the experimental — to rapidly fall to zero above 
the Debye energy. Very often this is not the case^i^i^ and 
a considerable background is found, which has been at- 
tributed to the presence of non-equilibrium phonons. It 
is however possible to correct for the background and to 
determine the appF(ey) function^. This method has 
allowed extracting the electron-phonon spectral function 
in many normal metals^, but can be applied also to su- 
perconductors above the critical temperature or driven 
normal by means of a magnetic field. Some examples 
will be discussed in Sect. IVII B 31 for the case of MgB2 
and in Sect. IVII Dl for the case of borocarbides. 



2. Thermal regime 

As opposed to the ballistic regime is the thermal (or 
Maxwell) one in which ^ <C a [see figure [5] (c)]. Some 
authors^ prefer to identify this regime by the condition 

^ a to make it explicit that electrons can undergo in- 
elastic scattering in the contact region as they normally 
do in the bulk. In this case, the resistance of the junction 
(already calculated by Maxwell) depends on the resistiv- 
ity of the metal^: 
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Joule heating occurs in the contact region and causes a 
local increase in temperature. The maximum tempera- 
ture Trnax at the center of the contact can be estimated 
by using the following expression^ 



(8) 



where Tj^ath is the bath temperature and L is the Lorenz 
number. In this case, at any finite bias the contact re- 
sistance is related to the resistivity of the material at 
Tmax > Tbath- Siucc in metals p increases with temper- 
ature, the I — V curves become S-shaped and the con- 
ductance decreases with bias^. Any spectroscopic infor- 
mation on the electron inelastic scattering is lost. Since 
the standard transport theory for bulk materials applies 
also to the contact, the FS is only slightly shifted in the 
direction of the electric field, as in Fig. [H^c). 



3. Intermediate regime 

Between the two aforementioned extreme regimes, the 
resistance of the contact can be expressed by a simple 
interpolation formula derived by Wexler— : 



R = 



2h 



(9) 



Here the first term is the Sharvin resistance and the sec- 
ond is the Maxwell resistance, multiplied by a function 
of the Knudsen ratio K. F is always of the order of unity. 
If the two metals are different (i.e. for a heterocontact), 
the resistance of the contact can be written a,SF^^ 
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assuming a spherical Fermi surface for both metals 1 and 
2. Here kp^min = niiii[/c_F,i, fcF,2] an 
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T = 



(vf,i + vpaY 



(11) 



In both eqs. [5] and [TO] the prevalence of the Sharvin or 
Maxwell term depends only on the size of the contact. 
For a junction between given materials, the Maxwell con- 
tribution dominates in large contacts, while the Sharvin 
one becomes more and more important on decreasing a. 

Between the thermal and ballistic regime one can also 
define the so-called diffusive regime in which the elastic 
mean free path £e of the electrons is small compared with 
the contact radius a but the diffusion length A = \/ lile 
for inelastic scattering is still bigger than a (a ^ A). The 
quasiparticles can now experience elastic scattering pro- 
cesses inside the contact region but not inelastic ones, as 
shown in the left panel of figure [IJb). The elastic scat- 
tering redistributes the quasiparticles isotropically over 
the FS, in an energy shell of width eV (right panel of 



Fig. mjb)). Though energy-resolved information is still 
available, the effective volume ^eff in which the inelastic 
scattering of electrons gives rise to the backflow current is 
now reduced by a factor of the order of a/ £ with respect 
to the ballistic regime (see eq. [S]). This is due to the 
fact that the probability for an electron to cross the con- 
tact, undergo an inelastic scattering event and then ffow 
back through the orifice is reduced by elastic scattering 
in the contact region. The intensity of the spectroscopic 
signal (proportional to —cPl/dV^) is thus strongly re- 
duced. Moreover, a different efficiency function must be 
used in the spectral function S{E) (see eq. S]), since the 
elastic scattering relaxes the requirement of momentum 
conservation. 



A. Determination of the conduction regime of a 
real point contact 

The radius of a real point contact (for example made 
by pressing a metallic tip against the sample surface) is 
unknown and, in general, experimentally inaccessible. As 
a matter of fact, the size of the actual contact is not re- 
lated to the apparent contact area or to the footprint of 
the tip® . So the problem arises of how to check whether 
the contact is ballistic or not. One possibility is to admit 
that the resistance of the contact Rn can be written as 
in the Sharvin formula, i.e. Rn = (4/0^)/(37ra^) where 
the product pi refers to the bank with the smaller Fermi 
energy (see eallOp and thus, generally, to the supercon- 
ductor. The condition a <ti £ can then be turned in a 
condition on the contact resistance: 
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(12) 



Alternatively, one can (very crudely) evaluate the contact 
radius a by means of 
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and then compare it to £. This estimation is based on the 
assumption that only one contact is present. In almost 
all real cases, because of the rather likely formation of 
parallel contacts, the value of a obtained in this way is 
nothing but an upper limit to the size of the contacts 
(whose number is unknown). As a matter of fact, in this 
case R]y is the resistance of the parallel as a whole and 
the resistance of individual contacts is necessarily larger 
than that. This means that, if a estimated from eq. [TO] 
is smaller than £, the contact (either single or multiple) 
is necessarily ballistic. If instead a > £, this does not 
necessarily mean that the contact is not ballistic. In these 
cases, the conductance curves [dl/dV vs. V) can help 
understanding what is the regime of conduction. If the 
conductance shows a downward curvature, for example, 
heating may occur in the contact. If the conductance 
shifts on increasing temperature, this may mean that a 
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Maxwell term (proportional to the resistivity) is playing 
a role. 

In the case of point contacts on superconductors, as 
we will see later on, some specific features show up in 
the conductance curves when the contact is not ballistic 
(see sect I VI A| . Moreover, a critical temperature of the 
junction smaller than the bulk can be due to a surface 
degraded layer but also, more banally, to Joule heating in 
the contact (so that the actual temperature of the contact 
Tmax is higher than that of the bath). 

IV. POINT-CONTACT 
ANDREEV-REFLECTION SPECTROSCOPY 
(PCAR) IN THE SUPERCONDUCTING STATE 

A. Andreev reflection 

Let us consider a normal metal (N) brought in direct 
contact with a superconductor (S), with no potential bar- 
rier between them. Let's apply to this junction a voltage 
V < A/e being A the energy gap in the S side. If the 
contact is ballistic, the whole voltage drop occurs at the 
interface. An electron coming from the N side will not 
be able to propagate through the interface because only 
Cooper pairs exist in this energy range in S. But if a 
hole is reflected and two electrons are transmitted in S 
as a Cooper pair (Fig. [3]) the total charge and momen- 
tum are conserved. This phenomenon is called Andreev 
reflections^ and can be theoretically described by solving 
the Bogoliubov-de Gennes equationsi^ at a N/S inter- 
face. The reflected hole has opposite wave vector and (if 
the Cooper pairs are singlets, as in all the cases analyzed 
here) opposite spin with respect to the incoming electron, 
so it traces back the trajectory of the incoming electron 
until a scattering event occurs. 

If the applied voltage is much greater than the gap 
{eV ^ A), all the electrons whose energy is lower that 
the gap still undergo Andreev reflection, but now their 
contribution to the current is constant and does no longer 
depend on the applied voltage. Instead, the electrons 
with energy higher than the gap are transmitted through 
the interface (see fig. ^ giving a voltage-dependent cur- 
rent. The total current for eV ^ A is thus^ 

V A 

I (xevpieV -A) + 2evF^~ h . (14) 

Rs eRs 

The second term of the right-hand side of eq. [T3]is called 
"excess current" and is the hallmark of the superconduct- 
ing state even at energies much higher than the gap. This 
result is exact only if the gap rises from zero up to the 
bulk value over a distance larger than the superconduct- 
ing coherence length ^. If the gap is instead modeled as 
a sharp barrier at the interface an additional term equal 
to A/3ei?s must be included. 

Because of Andreev reflection, the conductance of the 
junction turns out to be doubled for V < A/e. This 
clearly suggests a simple way to determine the energy 




FIG. 3: Electrical transport at an ideal (barrierless) N/S in- 
terface at T = 0. Incoming electrons with eV < A are re- 
flected as holes and, for each electron, a Cooper pair is trans- 
mitted (Andreev reflection) . Electrons approaching the inter- 
face with eV > A are normally transmitted as electron-like 
quasiparticles. 



gap in the S side by point contact spectroscopy. This 
technique is often referred to as point-contact Andreev- 
reflection spectroscopy (PCAR). From the solution of the 
Bogoliubov-de Gennes equations near a N/S interfaceii 
it is possible to note that Andreev reflection does not 
occur abruptly at the interface but over a length scale of 
the order of ^. In general ^ is also the length over which 
A is depressed due to the proximity effect generated by 
N on S. However, if the contact size is smaller than ^ this 
effect can be neglected. 

As already mentioned, PCAR requires that the gain 
in energy of the electrons crossing the junction is well 
defined. This is true in the ballistic regime but, also, 
in the diffusive one. If one wants to measure the gap 
by PCAR, it is clear that the voltage across the junc- 
tion will reach values of the order of, and even greater 
than, the gap A. If the contact is ballistic, using the 
value for the carrier density in the free-electron model, it 
is possible to show^ that the velocity of electrons across 
the contact, at F ~ A/e, is on the order of the depair- 
ing velocity in the superconductor. In other words, the 
current density becomes overcritical in the contact. Just 
outside the contact the current spreads out, its density 
decreases and will reach the critical value a short distance 
away from the actual junctio n^^i^^ , as shown in figure [3] 
(a). If the size of the overcritical region is smaller than 
the coherence length ^ the spectroscopy is still possible^, 
because superconductivity cannot be quenched over dis- 
tances smaller than ^. Therefore it is necessary to adopt 
contacts (see figure|4](a)) which are smaller than the elec- 
tron mean free path (to avoid heating effects) and smaller 
than the coherence length (to avoid proximity effect and 
destruction of superconductivity in the contact region)^. 
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■ Z = 0.1 

■ Z = 0,2 
• Z = 5 



FIG. 4: (a) Sketch of a ballistic Andreev-reflection point con- 
tact whose radius a is much smaller that the electron mean 
free path I and the coherence length ^; (b) A polar plot of 
the normalized transparency of an NS junction as function 
of the angle of injection of the current for different values of 
the Z parameter. 



The Blonder, Thinkam and Klapwijk (BTK) 
model 




2 2 

eV/A 



FIG. 5: (a) Normalized conductance curves a{E)/TN of a N/S 
interface at T = calculated within the BTK model (ea l20p 
as a function of the barrier parameter Z, from pure Andreev 
{Z = 0) to pure tunneling {Z — 10) regimes, (b) Effect of the 
thermal smearing on the normalized conductance (note that 
we used A = 3 meV and Z = 0.2 in all the curves). 



Even if Andreev reflection was discovered in the early 
60s, it was only in 1982 that Blonder, Thinkam and 
Klapwijk^*' (from now on referred to as BTK) gave a 
complete, even though simplified, theoretical discussion 
of the phenomenon, including the effect of a finite trans- 
parency of the interface. The most noticeable simplifica- 
tion is that the model is ID, i.e. all the involved momenta 
are normal to the interface and parallel to the x axis. The 
barrier is represented by a repulsive potential UoS{x) lo- 
cated at the interface, which enters in the calculations 
through the dimcnsionless parameter 



Of course, the smaller is Z, the more transparent is the 
barrier. The parameter Z is originally meant to represent 
the effect of the typical oxide layer in a point contact, the 
localized disorder in the neck of a short microbridge or 
the intentional oxide barrier in a tunnel junction. Ac- 
cording to the BTK model, calculated at T = 0, the elec- 
tron coming from the N side can undergo four processes 
whose probabilities are: 

A probability of Andreev reflection. The probabil- 
ity decreases with increasing Z for eV < A and is 
always small ioi eV > A; 

B probability of normal specular reflection. This 
probability increases with Z , i.e. on decreasing the 
barrier transparency; 

C =^ probability of transmission in S as electronlike 
quasiparticle (ELQ). The probability decreases if 
Z increases but it is always zero for eV < A; 

D ^ probability of transmission with FS crossing (i.e. 
as holelike quasiparticle, HLQ). The probability is 
small for eV > A and always zero for eV < A. 



Of course the sum of the four probabilities must be equal 
to 1. Fig [3] shows the particular case of a barrier less {Z — 
0) N/S junction at T = 0, where only the terms A and 
C are present. 

In can be shown that the expression of the total current 
across the junction, at T=0, is given by^ 

/oo 
[f{E - eV) - f{E)][l + A{E) - B{E)]dE 
-OO 

(16) 

where f{E) is the Fermi distribution function, A{E) 
and B{E) are the coefficients giving the probability 
of Andreev and ordinary reflection, and the quantity 
[1 + A{E) — B{E)] (which is the transmission probability) 
is often indicated by cr{E). Note that, although (t{E) 
is formally written only as a function of A and B, the 
contribution of C an Z? has been taken into account in 
the calculations. Iq is a constant which depends on the 
area of the junction, on the density of states and on the 
Fermi velocity. The derivative of the current with re- 
spect to the bias, dI]\fs/dV, provides the conductance of 
the junction. When divided by the conductance of the 
same junction when the superconductor is in the normal 
state, dl^N/dV, this gives the normalized conductance 
of the junction, G (which is the outcome of PCAR ex- 
periments). 

Here, instead of giving the explicit expressions for the 
probabilities A, B, C and D that can be found easily in 
literatur e'^'^i^" , we prefer to show in detail the results of a 
different approach^! that allows writing the AR normal- 
ized conductance at T = 0, G = {dl / dV)NS / {dl / dV)MN 
as a function of the quantities Nq{E) = E/y/E"^ — A^ 
and Np{E) = A/^/E'^ - A^ whose real parts arc the BCS 
quasiparticle and pair density of states, respectively. 

We can start from the deflnition of the transparency 
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T^v of the barrier in the BTK approximation of current 
injection totally perpendicular to the N/S interface: 



1 



1 + Z2 

and then we introduce the function: 



I E - y/E^ - A2 E - - A2 



E + ^/W^^^ 



It is trivial to show that: 



l{E) 



N,{E) - 1 



(17) 



(18) 



(19) 



Note that 7(i?) is a complex function even if the gap 
A is real, as in the BCS case, since Np{E) and Nq{E) 
become imaginary for E < A. By using these definitions 
it is possible to demonstrate that the BTK conductance 
at T = is given by: 



aiE) = TN 



l + rNh{E)\^ + {rN-l)hiE)^\' 
|l + (rjv-l)7(i?)2p 



(20) 



The calculated normalized conductance G{E) = 
a{E)/T]\f is shown in Fig. [HKa) for various values of Z 
and for A = 3 meV. In a perfectly transparent junction 
(Z = 0, pure Andreev regime) the conductance within 
the gap {\eV\ < A) is doubled with respect to the normal- 
state one. When Z > 0, two peaks appear at \eV\ « A 
and their amplitude increases on increasing Z while the 
zero-bias conductance (ZBC) is depressed. Finally, at 
Z > 10, the normalized conductance at T = coincides 
with the BCS quasiparticle density of states, i.e. the real 
part of Nq{E). Indeed, it can be demonstrated that the 
results of the BTK model for Z ^ oo coincide with the 
standard results of the theory for NIS (I=insulator) tun- 
nel junctions. Hence, the BTK model can reproduce, by 
simply changing a parameter, all the different experimen- 
tal situations corresponding to different transparencies at 
the N/S interface, from zero to infinity. 

Equation [20] is particularly useful to discuss the exten- 
sions of the simple BTK formalism we will present in the 
following sections. As a matter of fact it should be borne 
in mind that, even if widely used as a simple tool for fit- 
ting the experimental PCAR spectra, the original BTK 
model is based on a large number of approximations and 
simplifications, i.e.: 

(1) All the calculations are made at T = 0; 

(2) The problem is ID, i.e. the current injection is only 
perpendicular to the plane interface; 

(3) The barrier is ideal and presents a null thickness; 

(4) The Fermi surfaces of both materials in N and S 
sides are spherical; 

(5) The Fermi velocities are the same in both sides; 



(6) The superconductor is supposed homogeneous and 
isotropic. Because of the mono-dimensionality, the 
gap A entering the equations is actually the gap in 
one single direction and represents "the" gap only 
if the order parameter is isotropic (i.e. it has a 
s— wave symmetry). 

(7) The N/S interface is atomically flat (somehow im- 
plicit in the ID current injection). 

In the following we will show that most of these restric- 
tions can be easily relaxed giving a more realistic tool for 
the analysis of PCAR experiments in a variety of uncon- 
ventional superconductors. 



V. BEYOND THE BTK MODEL 

1. Finite temperature 

The calculation of the differential conductance of a N/S 
junction at finite temperature is a quite easy task. It can 
be simply accomplished by introducing in the equation 
for the current the standard convolution with the Fermi 
function at finite T, f{E, T), and then taking the deriva- 
tive of the current with respect to the bias voltage, i.e.: 



dl, 



NS 



dV 



iV) 



d f 

lo^ J [fiE-eV,T)- f{E,T)]a{E)dE 



(21) 

where a{E) is given by eq. [201 In figure [Sfb) the ef- 
fect of the thermal broadening on the normalized conduc- 
tance is calculated by using a temperature-independent 
gap A = 3 meV and Z = 0.2. At the increase of T the 
two peaks typical of the AR at Z ^ are smeared out 
finally leaving a single zero-bias maximum at ^ > 10 
meV. If the (supposed BCS) temperature dependence of 
the gap A(T) is taken into account in the expressions of 
Nq{E) and Nc{E), i.e. in the (j{E), the curves become as 
shown in Fig[6| (b) . The AR features now correctly dis- 
appear at the critical temperature of the contact (usually 
equal or very close to the Tc of the superconductor) . 

The pre-factor Iq of eq[2T] is expressed in terms of the 
normal density of states of the two materials and thus 
could, at least in principle, depend on temperature and 
on energy: in this case it should be brought inside the 
integral, and would no longer simplify when normalizing. 
This could be the case when the normal state conduc- 
tance is found experimentally to change with tempera- 
ture or to be voltage-dependent, as it is in cuprates^ and 
in the recently discovered Fe-based superconductors (see 
sect I VII C|) . However, one usually assumes for simplic- 
ity that Iq is constant and uses the expression for the 
normalized conductance to fit the experimental PCAR 
spectra. From the experimental point of view, however, 
these cases present the extra problem of defining what is 
the normal-state conductance to be used for the normal- 
ization, as we will show in sect lVII Cl 
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2. 2D or 3D BTK model 

If the current injection was really only perpendicular 
to the interface as the BTK model assumes, one could 
in principle probe the k dependence of the gap by mak- 
ing directional PCAR (DPCAR) measurements on the 
different crystallographic planes of high-quality super- 
conducting single crystals. Actually, charge carriers can 
approach the interface from any direction and the only 
condition set by the AR theory is that the component of 
the k vector parallel to the interface is conserved in all 
processes. This implies, for example, that the reflected 
hole comes back in N with k opposite to that of the in- 
cident electron and traces back its trajectory until the 
first scattering event in N occurs (see Fig. |l{a)). In the 
S side a Cooper pair propagates essentially in the same 
direction as the incident electron (neglecting the small re- 
fraction due to the expansion of the FS). Calling 6* at the 
angle between the direction of the incident electron and 
the normal to the interface, the conservation of trans- 
verse momenta leads to the following dependence of the 
transparency rjv on 6^: 



tn{On) 



cos{9n)'^ 



(22) 



Of course eql22l coincides with eq[T7]for On = Q- In fig- 
ure IHb) the angular dependence of the normalized trans- 
parency (i.e. tn{6n)/tn{0)) is shown for different values 
of Z. When Z = all the quasiparticles are transmitted 
with the same unitary probability in the whole half-space 
— 7r/2 < 9n < 7r/2, but at the increase of Z the transmis- 
sion becomes progressively weaker and more directional 
around the perpendicular to the interface. Strictly speak- 
ing, the injection is always in the whole half-space but 
one can decide to conventionally fix a threshold (e.g. 75 
% of the maximum transparency) to determine an equiv- 
alent injection angle 9*. In the limit Z > 10 (tunnel 
regime) one gets 9* « ±30°, i.e. the tunneling process 
is certainly highly directional. For the typical Z val- 
ues observed in real PCAR experiments (~ 0.2 — 0.5), 
6* ranges between ±70° and ±52° thus evidencing the 
reduced directionality of the PCAR technique. In ad- 
dition to these "theoretical" limitations, some practical 
problems have to be taken into account. Irrespective of 
the way the PC are realized (needle-anvil or "soft" tech- 
nique) the contact footprint has a relatively large area 
(some hundreds of square microns). If this area con- 
tains crystal-growth terraces, defects, pits or cracks the 
probability to have some contacts along a different crys- 
tallographic direction becomes high. Directional PCAR 
(DPCAR) spectroscopy can give reliable results only if 
very high-quality single crystals with highly regular (and 
large) surfaces parallel to the crystallographic planes are 
used. Despite these limitations, we will show in the ex- 
perimental survey (sect. IVII[) that recent DPCAR exper- 
iments were able to precisely determine the anisotropic 
properties of the gap in several unconventional supercon- 
ductors. 



As shown in eq. [22] the barrier transparency depends 
on the direction of the incoming electron in the N side. 
By introducing this expression in eq |201 integrating over 
the whole half-plane and properly normalizing, we get 
the normalized conductance at T = 02i: 



G2d{E) 



J^^ cr{E,9N)cos9Nd9_ 



N 



J^i tn {9 n) COS On d9i\f 



(23) 



The calculation of G2D at any temperature can be done 
as in ea l211 by a convolution with the Fermi function. 

When the system has rotational symmetry around the 
axis normal to the interface (i.e. the gap is isotropic 
and the FS is spherical) this approach can be consid- 
ered as the 3D extension of the BTK model. Figure [Sl^a) 
shows the comparison of two normalized conductances 
at r = and T = 4 K calculated with the standard ID 
BTK model and with its 3D version. The angular inte- 
gration leads to a remarkable depression of the AR signal 
when < Z < 10. Obviously, when Z = (completely 
transparent junction) or Z > 10 (tunneling regime) the 
two approaches yield the same results. In figure Elb) a 
complete temperature dependency of the normalized con- 
ductance calculated by using the 3D model and assuming 
a BCS A(T) dependence is reported. It is trivial to show 
that the 3D normalized conductance practically coincides 
with the ID one calculated for a properly enhanced Z 
value. Probably this fact explain why the standard ID 
model is still largely used in fitting the experimental data. 
Nevertheless, problems can arise when comparing the Z 
values obtained by the two different approaches, partic- 
ularly in the cases where the value of Z has remarkable 
consequences on the interpretation of the physical pro- 
cess occurring at the interface, as, for example, in the 
study of ferromagnet-superconductor PCAR junctions. 



3. Fermi velocity mismatch at the interface 

In a realistic system the Fermi velocities will be dif- 
ferent on the two sides of the contact. The mismatch 
of the Fermi velocities gives rise to carrier reflections at 
the interface even when no barrier is present. This effect 
was initially introduced in the original BTK theory^ by 
adopting an effective barrier parameter: 



-'eff 



Z2 



4r 



(24) 



where r = vpi/vp2 is the ratio of the Fermi velocities in 
the superconducting and in the normal side. The normal- 
state resistance at high voltage is given hy Rn = i?s(l ± 
Z'^j^f) where Rs is the Sharvin resistance®. 

In the 3D version of the model^i the situation is more 
complex. To account for the possibility of different ef- 
fective masses in N and S, the parameter r of eql21| is 
replaced by Aq = ks/k^r. The "refraction" of quasi- 
particles at the interface is due to the conservation of 
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FIG. 6: (a) Normalized conductance curves calculated at T = 
and T = 4 within the ID BTK model^S (dashed lines) and 
within its 3D generalization^i (solid lines) using A = 3 meV 
and Z — 0.2. (b) Temperature dependence of the conductance 
curves calculated within the 3D BTK model with Z — 0.2 and 
assuming for the gap a BCS temperature dependence with 
A(r = 0) = 3 meV, n = 19.73 K. 



transverse momentum, i.e. sin{9N) = ^osin{9s) where 
On and 9s are the incidence and transmission angles, re- 
spectively. Under these conditions it is possible to show^i 
that the normal transmission probabihty (eq. I22[) be- 
comes: 



tn{9n,9s) 



4Ao cos 9n cos 9s 
[cos Sat -I- Ao cos 



(25) 



By introducing this expression in the formula for the 
superconducting transmission probability (eq |20p and ex- 
pressing ^5 as a function of by using the "refraction" 
relation sinifiN) = ^asin{9s), one formally obtains the 
same expression for the normalized conductance G2d{E) 
as in eq. [221 that now, however, accounts for the mis- 
match in the Fermi velocities. Incidentally, when Ao < 1, 
i.e. kjv > ks, a "total reflection" of electrons occurs at 
the interface for injection angles \9]\[\ > sin~^XQ. In this 
case the integral in bas to be restricted to this limit 
angle^l. It seems that the condition Aq < 1 could eas- 
ily apply in the case of a superconductor with a small or 
very small FS and, thus, this problem could be important 
in new unconventional superconductors. In the opposite 
case, Aq < 1, can vary in the whole half-plane while 
the range of Os is restricted. Anyway, whatever the ap- 
proach to the problem is, it turns out that the global 
effect of a mismatch of Fermi velocities at the interface is 
simply described by a sort of "renormalization" of the Z 
values of the kind described in eq. [M] As a consequence, 
apart from extreme and hypothetical cases showing very 
large (or very small) Aq values, the effect of the mismatch 
cannot be separated from the standard experimental vari- 
ability of Z values, unless one is able to determine the 
true Z value at the interface. 



4- The broadening parameter 

Even if the BTK model allows a correct inter- 
pretation of some experiments in low-temperature 
superconductors^, in most cases it predicts much sharper 
gap features than those actually observed in the low- 
temperature conductance curves. This means that the 
AR structures in the experimental spectra are not only 
depressed in amplitude but also spread in energy. This 
effect can be attributed to the reduction of the quasi- 
particle lifetime, resulting from: i) the imaginary part of 
the quasiparticle self-energy. This term is "intrinsic" but 
very small, as discussed in the tunnel regime by Dynes et 
al-^ ; ii) inelastic quasiparticle scattering processes occur- 
ring near the N/S interface (surface degradation, contam- 
ination etc. either at the N or the S side)S^. This term 
is "extrinsic" and much larger than the previous one. By 
properly solving the Bogoliubov-de Gennes equations in 
the presence of an inelastic scattering term, it has been 
shownii^ that it is possible to globally take these ef- 
fects into account by including into the BTK model a 
single broadening parameter F in the form of an imagi- 
nary part of the energy, i.e. E E + iT. F can thus be 
considered as the sum of the "intrinsic" lifetime param- 
eter Ti = h/Ti and the "extrinsic" one Fg = H/Te, being 
Ti^e the corresponding intrinsic and extrinsic lifetimes. 
There is actually a third possible origin of broadening 
of the conductance curves that can be accounted for by 
using F, i.e. a distribution of gap values (in anisotropic 
superconductors). In this case, F simulates the effect of a 
convolution of the theoretical conductance with the gap 
distribution (an example is presented in Sect I VII D|) . 

Introducing F in the BCS quasiparticle density of 
states leads to the modified expression^SiS^: 



N{E,T)^^ 



E + iV 



A2 



(26) 



F enters the BTK model or its generalizations through 
Nq{E) and Np{E) in eq. [HI thus modifying a{E,9M) 
and the conductance G2d{E) (eq. [23|l . Fig. [Tj^a) de- 
picts the normalized conductance G2d{E) calculated us- 
ing Z — 0.25 and different values of the ratio F/A. 
The broadening effect of F cannot be reproduced by any 
combination of parameters of the standard BTK theory 
unless one convolutes the zero-temperature conductance 
with the Fermi function at a fictitious temperature higher 
than the actual one. This approach is sometimes implic- 
itly used indeed when the experimental smearing of the 
curves is treated in terms of a Gaussian broadening. Such 
a procedure is not theoretically founded and mixes the 
actual thermal smearing with the other broadening ef- 
fects, which are instead well distinct. Finally, even if it is 
common (and reasonable) opinion that the best conduc- 
tance curves should allow a fit with F/A < 0.5, large F 
values might be sometimes necessary (for example in the 
presence of a wide gap distribution). This does not nec- 
essarily prevent the determination of the gap by means 
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of a fitting procedure, which is indeed possible even when 
r/A > 1 (especially if Z is sufficiently large). 



5. Energy dependence of the order parameter 

It is well known that the mean-field BCS definition of 
a constant superconducting order parameter A is only 
a crude approximation of the physical reality. Actu- 
ally, even in the weak-coupling regime A is a func- 
tion of the energy and shows a small energy-dependent 
imaginary part. The signatures of this energy depen- 
dence on the normalized tunneling (or AR) conductance 
curves are extremely small but, when the intensity of 
the electron-phonon coupling increases (strong-coupling 
regime) they become visible. By solving the Eliashberg 
equations for the strong-coupling regime starting from 
the electron-phonon spectral function a^F{E) and the 
Coulomb pseudopotential /i* (direct solution) it is pos- 
sible to obtain the full energy dependence of the order 
parameter 1\{E) = 'St/\{E) + i'^A{E). The imaginary 
part of A{E) increases at the increase of the coupling 
and accounts for the finite lifetime of Cooper pairs. By 
introducing the function A{E) into the expression for the 
quasiparticle density of states (eq. [26] with F = 0), small 
deviations from the BCS DOS at the typical phonon 
energies are observed, due to the electron-phonon in- 
teraction (EPI). It is well known that also the inverse 
procedure works (but only approximately in multi-band 
superconductors!^'*) i.e. starting from the EPI structures 
in the experimental tunneling conductance it is possible 
to obtain a^F{E) and ^* by the inverse solution of the 
Ehashberg equations. 

Since the BTK theory (and its modifications discussed 
so far) coincides with the BCS theory for superconduct- 
ing tunnel in the limit of large , it is easy to predict 
that the introduction of A(£^) into the BTK expressions 
will lead to EPI structures in the normalized conductance 
for any Z value in the ballistic regime. This is indeed the 
case, as it can be explicitly demonstrated'^^. A simpli- 
fied approach to the problem was presented in Ref.^^, 
where simple asymptotic expressions for the normalized 
conductance at eV ^ A in the tunnel {Z oo), ballistic 
and diffusive regime were obtained by taking into account 
phonon self-energy effects on the order parameter. Let 
us instead show here an example of the complete pro- 
cedure applied to a "classic" strong-coupling supercon- 
ductor. First, we calculated the A{E) function of lead 
starting from its EPI spectral function (top curve in Fig. 
Wl,c)) and assuming ^* — 0.11. A{E) was thus introduced 
in the expressions of Nq{E) and Np{E) finally leading to 
the point-contact normalized conductance shown in Fig. 
[Tljb) for different Z values. As expected, the normalized 
conductance at eV < Apb coincides with the standard 
BTK one^. At eV Apt + Epu (where Eph repre- 
sents the range of energies where a'^F{E) pb ^ 0) the EPI 
structures appear for any Z value but their amplitude in- 
creases with Z. Fig. shows the sign-changed first 
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FIG. 7: (a) Conductance curves calculated at T = IK within 
the 3D BTK model by using A = 3.0 meV, Z = 0.25 and 
increasing values of F. (b,c) Normalized conductance and its 
voltage derivative {—dG/dV) calculated for different Z values 
within the 3D BTK model, in the case of Pb, by using the 
energy-dependent order parameter A{E) as obtained from the 
electron-phonon spectral function a^F{E) (top curve in panel 
(c), here shifted in energy for ease of comparison). 



derivative of the normalized conductance —dG/dV = 
—d'^lNs/dV'^ vs. V compared to the a'^F{E)pb (top red 
curve). Even if the EPI structures shift to higher ener- 
gies and their amplitude is depressed at the decrease of 
Z, the use of DPCAR spectroscopy in very high-quality 
single crystals to access quantitative information on the 
a'^F{E) and its dependence on direction, temperature 
and applied magnetic fields proved to be a feasible task^^. 



6. Anisotropic order parameter 

The assumption of an isotropic (s-wave) order param- 
eter (OP) makes the BTK model particularly simple, but 
this constraint must be relaxed if one wants to describe 
systems in which the OP is instead anisotropic, i.e. it de- 
pends on the wavevector k in the reciprocal space. This 
happens for example in high- Tc cuprates, where at least 
one component of the OP has a d-wave symmetry^. Gen- 
erally speaking, the anisotropy of the OP can have two 
different origins: (i) the OP has a true k dependence (at 
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least along some planes of high symmetry) on the sin- 
gle FS sheet where it opens; (ii) different isotropic OPs 
open on different sheets of the FS of a multiband system. 
Strictly speaking, in this case the OP is not anisotropic 
but appears so when it is measured by techniques with 
null or poor resolution in the k space. Of course, more 
complex cases with multiple anisotropic gaps can in prin- 
ciple occur, which could be probably elucidated only by 
experimental techniques with full k-space resolution (e.g. 
high- resolution ARPES). In this section, we will show 
how to account for a single anisotropic OP within the 2D 
BTK model. The more complex effect of multiple OPs 
on different FS sheets and the influence of the shape of 
the FS itself will be addressed in the next section. 

The problem of introducing the OP anisotropy into the 
expression of the superconducting transmission probabil- 
ity a{E, 9n) was solved in Ref.— in the most general case. 
Here we will give a simplified "operative" description of 
the general results. Let us suppose for simplicity that the 
OP has a k dependence only in the kxky plane and that 
X is the direction normal to the flat junction interface. 
Let the system have a translational invariance along the 
kz axis so that the problem reduces to a two-dimensional 
one, i.e. the FS is a cylinder. We also suppose that the 
current injection occurs in the plane xy (a5-plane con- 
tact) and that Aq = 1, i.e. there is no refraction of quasi- 
particles at the interface and both the integration angles 
6n and 6*5 span in the range [— 7r/2, 7r/2]. Let the OP 
A be a function of the angle 63 with which electron-like 
quasiparticles (ELQ) are injected in S. The specific ex- 
pression of A{9s) depends on the kind of symmetry the 
OP shows in the k space. To take into account the pos- 
sible rotation of the crystallographic a axis with respect 
to the normal to the interface {x axis) we also introduce 
the angle a [see figs. [HJa) and[5{b)]. Since ELQ and 
HLQ are injected in S with angles 9s and —Os, respec- 
tively, they feel different OPs, namely A_|- = A{6s — a) 
(for ELQ) and A_ = A(-6's-a) (for HLQ). Under these 
conditions, the superconducting transmission probability 
becomes^i: 



g\ anisotropic s 



|1 + (tjv 



where 



(27) 



and — {ip_ — ip+), ip± being the phases of A±. When 
A± are real quantities, then their phase can only be ei- 
ther or TT and the same holds for (pd- The choice of 
a determines the 9s intervals in which the phase differ- 
ence (fid is or TT. If A± do not show sign changes as 
a function of 9s, then (fd = independently of a. tn 
appearing in eg 1271 has the same expression shown in eq. 
[25l Putting (j{E,9m) in eq. [23] one finally obtains the 
total (integrated) normalized conductance at T = 0. The 
convolution with the Fermi function as in ea l21l will fi- 
nally give the theoretical curves to be compared with the 
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FIG. 8: (a,b) Polar plot of the magnitude of the OP in the case 
of anisotropic s-wave symmetry (a) and d^2_y2 symmetry (b). 
The sign of the OP is indicated by -|- and — . The angle a 
(here equal to 7r/4) between the crystallographic a axis and 
the X axis (normal to the interface) is also shown . (c,d) 
Normalized conductance curves at T = IK calculated within 
the generalized 2D BTK model^i with a = for different 
values of Z and in the same OP symmetries as in (a, b). 



experimental results at any T. FigurelHJc) shows the nor- 
malized conductance at T = 1 K for different Z values 
and a = in the case of anisotropic s-wave symmetry of 
the pair potential, where A+ = Ai -I- A2 cos^[2(05 — a)], 
A_ = Ai-|-A2C0s4[2(-6i5-a)] (Ai = 1.5 meV, A2 = 1.5 
meV) and ifd = 0- Figure [UJd) shows the normal- 
ized conductance for the same values of the parameters 
in the case of dj;2_y2-weLve symmetry of the OP, where 
A+ = Ai cos 2(6*5 -a), A_ = Ai cos2(-6is-a) (Ai = 3 
meV). In both cases the shape of the normalized conduc- 
tance is quite different from the behavior shown in the 
s-wave case for the same Z values [figure [5] (a)]. In par- 
ticular: i) the anisotropic s curves show a four-peak (or 
two-peak and two-shoulder) structure similar to that ob- 
served in MgB2 ii) the 6,^2 _y2 curve with Z = b presents 
the well-known V-shaped conductance at low bias, while 
the one with Z = shows a cusp at zero bias instead of 
the flat region typical of the s-wave symmetry. In the d- 
wave symmetry, a = ±7r/4 gives Lpd = it for any value of 
^5 and the normalized tunneling conductance (for high 
Z) presents the well-known zero-bias conductance peak^. 
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7. True shape of the Fermi surfaces and momentum 
dependence of the pair potential 

Taking into account in the calculations for the PCAR 
conductance the true shape of the Fermi surfaces in N 
and in S, the possible k dependence of the pair poten- 
tial and the possible existence of multiple sheets of the 
FS - where the OP can assume different values - is a 
rather complicated task, from both the conceptual and 
the numerical point of view. Let us proceed step by step 
following the approach reported in^liS^. We will neglect 
possible interference effects between bands, that can lead 
to the formation of bound states at the surface as dis- 
cussed in Reff^. 

First of all the materials used in the N side are usually 
good conductors (Au, Ag, Pt, Cu, Al) for which the ap- 
proximation of a spherical FS is reasonable. So here we 
restrict the analysis to the shape of the FS in the super- 
conducting material. In the most general case the FS is 
divided into different sheets. Let us label them with the 
subscript i and call n the unitary vector in the direction 
of the total injected current, perpendicular to the con- 
tact interface. As a consequence the components along 
the direction n of the Fermi velocities at wave vector k in 
the ith FS sheet of the superconductor are v^k-n = fik.n 
where v^k — \ {Vk [i?i(k)]}. Of course, due to the pre- 
vious approximation, the corresponding quantity in the 
normal metal is viq.n — • n, being vjv the (constant in 
magnitude) Fermi velocity in the normal material. The 
ith component of the total current flowing through a per- 
fectly transparent {Z = 0) interface with no mismatch of 
the Fermi velocities (Aq = 1) in a ballistic PCAR exper- 
iment on a superconductor with isotropic OP is thus^: 



ik,n)FSi 



N.kV,lr,ndSl. 



(28) 



FS, 



where N^kiEp) = 1/ {47r3|Vk [Et(k)] |}^^ is the density 
of states of the ith band at the Fermi energy and wave 
vector k in S, dSp is the elementary area on the FS in 
S and QpSi is the integral over the ith FS sheet. The 
integral in eq. [28]is limited to values Wik.n > 0. Obviously 
Si^n has the meaning of area of the projection of the 
ith FS sheet along the n direction, i.e. on the interface 
plane perpendicular to n. It is the area of the ith FS 
sheet of the superconductor "seen" along the direction 
n. Of course, under these restrictive conditions every 
contribution to the total conductance from the ith FS 
sheet can be evaluated by using the same kind of integral, 
i.e. is proportional to the projected area Si n. It means 
that the total conductance "seen" along the direction n 

is (cr(£'))/||„ = Y.i<^i{E){^i'^'"iKn) FSi = Y.i^i{E)Si,n 

and the total normalized conductance is: 



{G{E)) 



/||n 



I],;(A^jkWik,n)_FS, 



where (Ji{E) is the BTK superconducting transmission 
probability (eq. [^0)1 of the ith FS sheet. In the case of 



different OPs on the different sheets of the FS the 
total normalized conductance will be dominated by the 
contribution of the (Ji(E) that corresponds the largest FS 
projected area along the n direction. As a consequence, 
directional PCAR experiments at Z = and Ao = 1 
can give information on the distribution and values of 
the isotropic OPs on the different FS sheets in a multi- 
band, multigap superconductor. It is quite obvious to 
expect a similar results also in the more general case of 
an anisotropic OP and of Z 7^ and Aq 7^ 1, but the 
calculation of the normalized conductance is now much 
more complex. 

First of all, if the OPs on the FS sheets are anisotropic, 
i.e. = Aj(k) = A^k (but stiU Z = and Aq = 1), then 
the superconducting transmission probability becomes a 
function of k and cannot be anymore extracted from the 
integral over the FS. The total normalized conductance 
thus becomes: 



(G(i?))/||n - 



Y.i ("'ik (-E) Af»kt'»k,«) FSi 
J2i{^ikVik,n)FS, 



(29) 



but us- 



where aiu_{E) is always expressed by eq. 
ing fu nctions N ^^^jE) = E/^E^ - A'f^ and Nf^{E) = 
Aik/V-B^ — A^^ that substitute the standard ones in the 
definition of j{E) in eq. [121 If the barrier has a finite 
transparency and there is a N/S Fermi velocity mismatch, 
the normal transmission probability of the barrier r/v is 
no longer identically 1. According to the standard 2D ex- 
tension of the BTK model shown before, tn (which here 
we call T for simplicity of notation) is given by eq. I25l that 
can be conveniently rewritten as a function of the pro- 
jections of the Fermi velocities along the n directional; 



T -ik ■ 



4Wik,nWAr,n 



(Wzk,n + WAr,„)2 -I- AZ'^vjf ' 



(30) 



By introducing this transmission probability inside the 
integrals over the FS both at numerator and denomina- 
tor of eq. [21] and taking into account that iVik'^ik.ri = 
Vik,n/vi]i wc finally obtain the total normalized conduc- 
tance at T = in the most general case: 



t'.k[(l',k.„-Hl'N,„)2-H4Z2„^] 



2ZT^)FS, 



(31) 

where a subscript n has been added to the expressions 
of o'ik(E) and Z just to include the possibility to have 
different Z values along the different crystallographic di- 
rections, a thing that is often observed in DPCAR ex- 
periments. In the case of large Z (tunneling regime) the 
weighting factor inside both the FS integrals of eq. [3T] 
reduces to A^ikW^j^ „ and the calculations are simplified. 
As previously, the presence of isotropic OPs on every FS 
sheet allows extracting ai{E) from the integrals. This 
is the approach recently followed by Brinkman et al. in 
Ref;^, where the total normalized conductance of MgB2 
has been written as a weighted sum of the partial con- 
ductances of the a and tt bands using the squares of the 
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FIG. 9: Left: Distribution of the pair potential values over 
the three different sheets of the FS of CaCg obtained by first- 
principle calculations?^ (by courtesy of A. Sanna and S. Mas- 
sidda). Right: theoretical normalized conductance at T = 
for current injection along the a axis {Za = 0.75) and along 
the c one [Zc = 1). The values of Z are taken from PCAR 
experiments. 
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plasma frequencies along the different crystallographic 
directions as weighting factors. Of course, independently 
of the isotropic or anisotropic properties of the OPs, if 
the current injection in a point-contact (or tunneling) ex- 
periment was a fully directional process the gap should 
not be seen along that directions where the FS has a null 
projected area. Actually, as we have seen in the previous 
sections, this is not the case, i.e. only a partial direction- 
ality is always present, which depends on the Z and Ao 
values. This explains why c-axis tunneling experiments 
on superconductors with a quasi-2D FS (cylinder parallel 
to kz) actually are able to measure the gap averaged over 
the ah plane. If the gap value A^k and the Fermi veloc- 
ity Vi]s_ are known at any k point of the «th FS sheet by 
first-principle calculations or by high-resolution ARPES 
experiments, then eq. [3T] allows the calculation of the 
PCAR normalized conductance at T = for a current 
injection along any crystallographic direction. An exam- 
ple of the results of this procedure^S is shown in Fig. [SI 
which shows the distribution of the pair potential values 
over the three different sheets of the FS of CaCg obtained 
by first-principle calculations'^-^ (left panel) and the theo- 
retical AR normalized conductance at T = for current 
injection along the a axis {Za = 0.75) and along the c 
one {Zc = 1) (right panel). The theoretical curves of 
Fig. 9, when properly broadened by P values close to 
the experimental ones, turned out to reproduce very well 
the experimental DPCAR results in CaCg^^-, as it will be 
shown in Sect I VII El 



FIG. 10: (a) Temperature evolution of the conductance curve 
of a Ag-paste contact on a Mgo.85 AI0.15B2 single crystal 
(normal-state resistance Rn = 13 fi) featuring clear dips. All 
the curves but the bottom ones are vertically offset for clar- 
ity, (b) Upper curves (shifted for clarity): the theoretical 
conductance curve obtained within the diffusive model with 
Z = (solid line) is well reproduced by that obtained within 
the BTK model with the same parameters but with Z = 0.59 
(dashed line). Lower curves: conductance curve calculated 
in the diffusive model with Z — 0.5 (solid line) and its BTK 
fit (dashed line) that requires Z = 1. In all cases, the other 
parameters of the models are always the same: T = 0.4 K, 
A = 1 meV and F = 0.2 meV. 



VI. NON-IDEAL EFFECTS IN THE CONTACT 
A. Dips 

The PCAR differential conductance often shows unex- 
pected sharp dips at voltage values larger than the su- 
perconducting gap, but sometimes very close to it, as 
shown in Fig. [TUT a). These dips are related to the su- 
perconducting properties of the S electrode since they 
never show up in NN junctions, but the BTK theory is 
unable to reproduce them. On increasing temperature, 
they generally shift to lower energies and generally affect 
the shape of the gap structures, as shown in FigfTOTa). 
For example, they can make a broad maximum centered 
at zero bias look as a sharp zero-bias conductance peak. 

It is commonly accepted that these dips indicate a 
non-ideal conduction through the contact. The detailed 
mechanism leading to their emergence was studied by 
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Sheet et ali^ who measured the evolution of the PCAR 
spectra of various N/S contacts (made with the needle- 
anvil technique) on progressively reducing the diameter 
of the point contact by withdrawing the tip in small steps, 
and found indeed that the dips are related to the regime 
of conduction through the junction, but also to the bias 
current. As a matter of fact, according to Wexler's for- 
mula, the point-contact resistance is generally the sum 
of a Sharvin and a Maxwell contribution, whose relative 
weight depends on a. If a is large, the dominant term is 
the Maxwell one, which contains the bulk resistivity of 
the two electrodes (eq. [TUl) . As long as the current flow- 
ing through the contact is small, the resistivity of the su- 
perconductor is zero; however, when the current reaches 
the critical value (/c) in the S side, a normal-state region 
can be created in S close to the junction, as discussed in 
Sect llV Al If this happens, the resistivity of the super- 
conductor starts playing a role and enters eq llOl giving 
a sharp increase in the voltage across the junction and a 
dip in the differential conductance. The same mechanism 
can be described as being due to the sudden disappear- 
ance of the excess current. Numerical simulations of the 
conductance, obtained by summing the I — V curves of 
a ballistic contact (given by BTK) to those of a typical 
bulk superconductor, give indeed results in good agree- 
ment with observations- 

An alternative explanation of the dips as being due to 
proximity effect was given in Reff^. The idea is that, if a 
proximity layer with depressed order parameter Aprox is 
present at the interface, Andreev reflection is limited to 
energies eV < Aprox, while quasiparticles can enter the 
S side only when eV > A^uik- This gives rise to dips in 
the conductance curves, at energies between Aprox and 
Abuik , which also necessarily shift to lower energies on in- 
creasing temperature because of the temperature depen- 
dence of the gaps. Refi^ also provides a model for the 
fit of the conductance curves that requires Ahuik, Aprox 
and Z as adjustable parameters and can be generalized 
to include a broadening term T. 

Very often, when analyzing conductance spectra with 
dips, a BTK fit is done ignoring the dips. However, even 
if this procedure introduces only a small error in the de- 
termination of the gap when the dips are small, it has 
been shown'^^ that a considerable overestimation of the 
gap can occur when they become large. 



B. Diffusivity in the contact 

In Sect nil 3l we mainly discussed the effects of a dif- 
fusive contacts in the case of a N/N junction. In N/S 
junctions, the diffusivity in the contact has been theo- 
retically addressed by Mazin et aZ.^-t?! and turns out to 
affect only the Z parameter. For instance, the conduc- 
tance of a diffusive junction with a given barrier parame- 
ter Z can be fitted with a ballistic (BTK) model with an 
effective Z* > Z. This is shown in figure [TUf b) where the 
conductances obtained within the diffusive model (solid 



lines) are compared with those calculated with the stan- 
dard BTK model (dashed lines). All the curves are cal- 
culated for r = 0.2 mcV, A 1 meV and T = 0.4 K. The 
upper curves show that the conductance in the diffusive 
model with Z = is well reproduced by the BTK model 
with Z = 0.59. Analogously, the lower curves indicate 
that when Z — 0.5 is introduced in the diffusive model, 
the obtained conductance corresponds reasonably to that 
obtained within the BTK model, but with Z — 1. This 
conclusion is also, and even more, true at higher temper- 
atures and for higher values of the lifetime broadening, 
i.e. when the curves are more smeared out. 



C. Inelastic scattering in the vicinity of the contact 

The inelastic scattering due to some layer with differ- 
ent composition at the N/S interface has been clearly sin- 
gled out experimentally in ref^?*?- where ballistic Andreev- 
reflection measurements were performed in Cu-Pb junc- 
tions with and without a very thin (« 2 nm) Pt layer 
in between. The PCAR curves of the Cu/Pt/Pb junc- 
tions were shown to be more broadened than those of 
the Cu-Pb contacts, and were well fitted by the BTK 
model by systematically using larger F values ~ though 
giving a good determination of the gap amplitude (note 
that, already in the original paper by Plecenik et ali^, F 
was introduced in the BTK model to take into account 
exactly these effects). 

Something similar is likely to happen in the "soft" 
point contacts, whose normalized conductance curves 
show a reduced amplitude and a larger broadening than 
those obtained with the conventional needle-anvil tech- 
nique. To identify the scattering layer in this case, we 
carefully measured the temperature dependence of the 
resistivity of the particular Ag paint used for the con- 
tacts. We found a residual resistivity at low temperature 
of 0.34 mfi cm (about 10^ times higher than that of pure 
Ag), and an enormously increased slope of p(T) at higher 
temperature. The former indicates a huge contribution 
of intergrain connectivity to the resistivity, and the latter 
a drastic reduction of the inelastic mean free path on the 
grain surface, which could well give rise to the observed 
broadening of the conductance curves. It must be said, 
however, that a contribution from a layer at the surface 
of the sample cannot be completely ruled out, and is in- 
stead proved by the fact that a similar broadening has 
been observed also in some PCAR spectra taken with the 
needle- anvil technique. This will be further discussed in 
the experimental survey (see sections I VII B II VII C|) . 

D. Spreading resistance 

For spectroscopic measurements to be reliable, elec- 
trons must not lose a significant energy while traveling 
through the electrodes. If at least one of the electrodes is 
highly resistive, a so-called spreading resistance Rsp must 
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be considered in series with the contact resistance, and 
this results in a shift of the conductance peaks to higher 
energies, leading to an overestimation of the gap^^. Ac- 
tually, a spreading resistance R^p is always present but 
usually plays a role only in measurements performed in 
thin films, while in bulk or highly conductive samples it 
is much smaller than the contact (junction) resistance 
and can thus be neglected. In the case of "soft" point 
contacts, one can wonder whether the Ag paste between 
the Au wire and the sample surface can give a significant 
contribution to Rgp- Actually, the resistance of the Ag- 
paste spot (approximately modeled as a cylinder with a 
diameter of 50 ^m) is as small as 0.086 51 even if a (largely 
overestimated) thickness of 50 /Ltm is assumed. This value 
is clearly negligible when compared to the contact resis- 
tance that is usually in the range 5 — lOOfi (depending 
on the material under study). 



VII. POINT-CONTACT SPECTROSCOPY IN 
MULTIBAND SUPERCONDUCTORS 

A. Two-band model for superconductivity 

The first theoretical study of multiband supercon- 
ductivity dates back to the late Fifties when Suhl, 
Matthias and Walkei— generalized the BCS theory 
to the simple case of a superconductor with two 
overlapping bands. The corresponding BCS Hamil- 
tonian contains two intraband terms of the kind 



»,kT »,-ki »,-k'i j.k'T 



and two interband terms 



of the kind ^j^j^' '=I,kT'^I-kiS,-k'iS,k'T ('^here ij 



1,2 is the band index). Vij is the (constant in the BCS 
approach) averaged pairing potential which results from 
boson emission and absorption by an i-j process, minus 
the corresponding shielded Coulomb interaction. In the 
absence of interband coupling (Vij = 0), the two bands 
would be completely independent, each featuring its own 
BCS gap and critical temperature. In the opposite case 
(only interband coupling, Vu = 0) the critical tempera- 
ture is the same, but there are still two gaps unless the 
partial density of states is the same in the two bands 
(TVi = A^2)- In general, through interband coupling the 
band with the higher superconducting Tc raises the crit- 
ical temperature of the weaker, or even induces super- 
conductivity in a nonsuperconducting band. The criti- 
cal temperature is defined as fesTc = LUkBOoe^^^^"" 
where Ag// is the effective coupling constant and is sim- 
ply the maximum eigenvalue of the matrix A^ = VijNj 
where Nj is the density of states at the Fermi energy (per 
spin) in the jth band. 

Figure [TT] shows the temperature dependence of the 
(normalized) gaps as a function of the normalized tem- 
perature in the BCS two-band model in different cases: 
i) bands completely decoupled {Vij = 0, solid lines). The 
Tc's of the bands depend on the relevant intraband cou- 
pling; ii) weakly coupled bands (dashed lines). While 
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FIG. 11: Temperature dependence of the gaps Ai and A2 in a 
two-band BCS model, calculated in the cases of: no intraband 
coupling (solid lines); weak intraband coupling (dotted lines); 
strong interband coupling (dash-dot lines). The intraband 
coupling constants are arbitrary; here we used those for MgB2 . 



Ai follows the same standard BCS temperature depen- 
dence (but with 2A/fc_BTc > 3.53), A2 features a high- 
temperature tail and closes at the same Tc as Ai; iii) 
strongly coupled bands (dash-dot lines). The small gap 
still deviates from a BCS-like behavior but smoothly de- 
creases on heating, to finally close rather quickly at Tc- 
The gap ratios 2A/kBTc for the two gaps are greater 
and smaller than the single-band BCS value 3.53, respec- 
tively. As we will show in the following experimental sur- 
vey, PCAR measurements in multiband superconductors 
have provided examples of all these three cases. 



B. Magnesium diboride 

After the publication of the theory for two-band super- 
conductivity, some of its consequences on various mea- 
surable quantities were calculated and possible marks of 
multiband superconductivity were found in conventional 
materials like Nb ^^'^° . In 1980 a clearer experimental 
evidence of multiband superconductivity was found in 
Nb-doped SrTiOgii by means of tunnel spectroscopy. 
Despite the fundamental importance of the result, the 
very low transition temperature of this compound (a few 
hundred mK) made its experimental investigation rather 
demanding and prevented its study from becoming very 
popular. The situation changed completely in 2001 when 
superconductivity below 39 K was discovered in MgB2, 
which remains up to now the most known and the most 
studied example of multiband superconductor. MgB2 
has a layered structure with graphite-like, honeycomb B 
layers intercalated by Mg planes with hexagonal close- 
packed structure^. Its electronic structure includes four 
a bands originating from sp^-hybrid B orbitals, and two 
TT bands due to the overlapping of the residual pz or- 
bitals. The Fermi surface is made up of nearly-2D cylin- 
ders around the F — ^ line (due to the a bands) and a 3D 
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tubular network related to the tt bands^^. Superconduc- 
tivity develops in the a bands below = 39 K mainly 
because of their coupling to the £'29 phonon modes^, and 
is induced in the tt bands through interband coupling. 

The key role of this two-band-system structure was 
soon witnessed by the failure of all the conventional 
single-band theories in describing the phenomenology of 
MgEy**^'^^. An effective two-band model was thus pro- 
posed, in which the four bands were grouped into two 
band systems (ct and tt). The anisotropic effective cou- 
pling constant for superconductivity Ae// — 1.01 actually 
indicates an intermediate coupling regime which is best 
described by the Eliashberg theory— i^i^. The calcula- 
tion of the gaps within a two-band Eliashberg model^^ 
gave A,^ = 7.1 meV and = 2.7 meV (see FiglT^. 
Similar values can be obtained within a BCS approaches,. 

An interesting feature of multiband superconductiv- 
ity in MgB2 is the role played by impurity scattering 
in the intraband and interband channels. According to 
Anderson's theorem^, it can be shown^ that, at least 
for small impurity concentrations, the intraband non- 
magnetic scattering has no effect on Tc and the gaps. 
The interband scattering on the other side, has a 'pair- 
breaking effect and thus decreases the critical tempera- 
ture Tc- According to the two-band model, in the limit 
of very strong interband scattering (dirty limit) a com- 
plete isotropization is asymptotically achieved, and the 
two gaps assume the same value so that one single gap is 
actually observed (dotted line in FigfT^). This is often 
referred to as "gap merging" . According to Eliashberg 
calculations in Refi^, at low temperature Adirty = 4.1 
meV with a corresponding reduced Tc = 25 K. For the 
sake of completeness, FigfT2k also shows the results of a 
fully-anisotropic Eliashberg calculation, based on the ac- 
tual momentum dependence of the electron-phonon cou- 
pling calculated ab-initic^^. This approach gives two 
distinct and non-overlapping distributions of gap values 
with average 6.8 meV and 1.8 meV. The differences from 
the two-band model arise from details in the calculations 
that are not worth discussing here. In any case, all calcu- 
lations show that the gap values on the two band systems 
are sufficiently different to be distinguishable also exper- 
imentally. 

According to the discussion of Sect IV 7i the shape of 
the FS (and in particular of the quasi-2D cr-band sheets) 
suggests a dependence of the PCAR or tunneling spectra 
on the direction of (main) current injection. Brinkman 
et al^ calculated the conductance curves of an ideal 
MgB2-I-N junction with various barrier transparencies 
within the Eliashberg theory. They expressed the nor- 
malized conductance G of the junction as the (weighted) 
sum of the BTK contributions of the two band systems: 
G = WttGtt + (1 — Wt)Gct^'*^^^. As expected, the weight 
Wa- = 1— depends on the direction of current injection. 
For / II c (and parallel to the axis of the nearly cylindri- 
cal CT-band sheets) is no more than 1% so that only 
the small gap A-^ should give detectable structures in the 
conductance curve. For / || a6, is maximum and equal 



to 33 %, so that four peaks corresponding to the small 
and the large gaps A^ and A^- are found in the conduc- 
tance curves^^. Note that the theoretical values of Wa 
are referred to ideal tunnelling current injection; slight 
differences are expected in PCAR experiments where the 
angle of effective current injection as defined in Sect IV 2| 
can be considerably larger. 
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FIG. 12: (a): Theoretical predictions for the gaps in MgB2 
within the Eliashberg two-band modelSS. The dotted line 
refer to the "dirty" limit (i.e. strong interband scattering). 
Similar results are obtained in the BCS approach. On the 
left, the gap distribution given by a fuUy-anisotropic Eliasherg 
calculation are shown^. (b) Experimental results of PCAR 
in MgB2, from Ref.— (circles),— (squares),^'* (up triangles),— 
(left triangles), (diamonds),— (stars). Lines indicate the 
predictions of the two-band model. 



1. Determination of ttie gaps in MgB2 

The earliest PCAR investigations carried out in MgB2 
polycrystals gave evidence of a single isotropic (s-wave) 
gap. Schmidt et alf'^ obtained A = 4.3 4.6 meV, 
while Kohen et alf'^ measured a gap A — 3.8 4.0 meV 
in higher-resistance contacts, while in a lower-resistance 
junction a smaller gap (A = 3 meV) was found, with 
reduced = 29 K. Laube et alM obtained an accumu- 
lation of gap values around 1.7 meV and 7 meV but never 
observed both of them in the same spectrum. Plecenik 
et alM studied the Andreev reflection curves of MgB2/N 
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junctions obtained in different ways, wfiose fit with tlie 
modified BTK model (with F — 0.8 meV) gave a gap 
A = 4.2 meV. A discontinuity in the temperature evo- 
lution of the gap suggested the existence of parallel con- 
tacts in clean and dirty regions of the sample, with a 
gap As = Att = 2.6 meV closing at Tc — 38 K and a gap 
Adirty = 4.0 meV closing at Tc = 22 K, respectively. The 
absence of A^- was probably due to a preferred c-axis cur- 
rent injection. The presence of a degraded layer on the 
sample surface suggested in^ and^ was confirmed by 
PCAR measurements performed with electrochemically 
sharpened tips of different hardness^^, which showed in- 
deed a decrease in the height of the conductance peaks 
(from 1.8 down to 1.25) and an increase in F from zero 
up to 1.2 meV on decreasing the pressure in the contact 
region from about 0.6 GPa down to 0.1 GPa. Spectra 
taken with the "soft" pressure-less technique had a height 
of only 1.15, and their fit with a single-gap BTK model 
gave F ~ 3 and A^irty = 4 — 5 meV with a reduced Tc- A 
reduced Tc was found also in Refi^, together with an in- 
crease in F and Z on decreasing the barrier transparency. 
All these results indicate an extrinsic contribution to F 
from inelastic carrier scattering in the barrier, not eas- 
ily accountable for in the theoretical model, and possibly 
due to a degraded or reconstructed layer covering the 
sample which can be broken by a tip but remains intact 
when the pressure is small or absent--. Indeed, it was 
shown experimentall}*25 that this effect can be simply 
accounted for by increasing the broadening parameter(s) 
in the modified BTK model. 

With the improvements in the sample quality, spec- 
tra with multiple gap features were readily obtained in 
films and polycrystal s^^^ , that allowed a fit by the two- 
band BTK model. In principle, the fitting function con- 
tains seven parameters: the two gap amplitudes A^ and 
Att, the broadening parameters Fo- and F^, two bar- 
rier parameters and Z^, plus the weight (so that 
Wa- = 1 — Wjr) for a total of 7 parameters. Some au- 
thors decided to use only one Z for both bands^i^ but, 
owing to the different Fermi velocities in the two bands, 
keeping Z^ and Z^^ as independent parameters is more 
general. Some authors also take = Ftt or even re- 
place them with a convolution of the T = conduc- 
tance with a Gaussian of width ojS?-. Others (including 
us) prefer instead to calculate the conductance at the 
correct temperature, and add F^ and F^r as imaginary 
parts of the energy in the BTK equation?-^ to account 
for all the sources of broadening discussed in Sect IV 4l 
Despite the number of free parameters, reliable values of 
the gaps can be obtained. This is certainly true for A^r 
which is quite strictly determined by the energy position 
of the relevant conductance peaks. The same holds for 
Act when the relevant peaks are observable - that means, 
for / 11 ab^. When the structures related to the large gap 
are only smooth shoulders (as in c-axis contacts or in ab- 
plane contacts at higher temperature), the uncertainty 
on Act increases. The evaluation of this uncertainty is 
not straightforward, because of the complex expression 



for the conductance in the two-band BTK model and the 
number of parameters. Indeed, an automated fitting pro- 
cedure is destined to fail, and one has to manually search 
for the parameters that allow minimizing the chi-square 
or the sum of squared residuals (SSR). Once the "best" fit 
is found, a range of parameters that give "acceptable" fits 
must be determined. This can be done by fixing a level 
of confidence for the chi-square or allowing a percent in- 
crease in the SSR. Then, the fit has to be repeated many 
times by changing all the free parameters so as to find 
the maximum variation of the gaps compatible with the 
fixed limits. Several fits made independently by different 
people normally ensure a good estimate of this range. 
Fortunately, some physical constraints limit the range of 
variability of some parameters. For example, w^r, Z„ and 
should not depend on either the temperature and the 
magnetic field; the intrinsic (lifetime) part of Fct and F^r 
can increase with temperature, but their usually much 
larger extrinsic part, related to the interface properties, 
should probably not. 

FigfTWb) reports the experimental results of various 
PCAR experiments in MgB2. In all cases apart from 
RefjS a two-band fit was used. All the data sets approx- 
imately agree with each other, apart from the early data 
by Bugoslavsky in thin films which show a reduced T^. 
The error bars are indicated only for some data sets and 
clearly increase on approaching because of the thermal 
smearing of the gap features. Because of the same effect, 
one may wonder whether the two gaps really close at the 
same temperature, since at high temperature the spectra 
show only a broad maximum and the two-band fit could 
be questioned. A conclusive answer to this issue and in 
favor of the two-gap model in MgB2 was found already in 
2001 by Szabo et alS^, who performed PCAR measure- 
ments in polycrystalline samples (squares in Fig |12r b)) 
obtaining gap values in very good agreement with the- 
oretical predictions. They found that the application of 
magnetic fields to the junctions resulted in a much faster 
suppression of the 7r-band features with respect to the 
cr-band ones. At high temperature or in c-axis contacts 
where no cr-band features are apparent, the disappear- 
ance of the dominant 7r-band structures allows unveiling 
the underlying cr-band contribution, with the emergence 
of two much well resolved maxima related to Act even at 
T = 30K. 

The synthesis of single crystals large enough to be used 
for PCAR allowed a step forward in the experimental in- 
vestigation of multiband superconductivity in MgB2, and 
in particular a study of the anisotropy of the spectra^^) 
by controlling the direction of (main) current injection. 
The soft-PCAR technique allowed us to make the con- 
tacts either on the flat surface of the crystals (c-axis con- 
tacts in the following, according to the nominal direction 
of current injection) or on their thin (50-100 /im) side 
(a6-plane contacts), which is very difficult by using a tip. 

Figure [T3] shows two examples of conductance spectra 
measured in a6-plane and c-axis contacts whose normal- 
state resistance is indicated in the labels. Note that for 
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FIG. 13: "Soft" PCAR spectra in MgB2 single crystals for 
different directions of current injection. The experimental, 
normalized conductance curves for a o6-plane contact (upper 
panel) and a c-axis contact (lower panel) are shown (symbols) 
and compared to the two-band BTK fit (solid lines). In the 
lower panel, the single-band BTK fit is also reported (dotted 
line). The normal-state resistance of the contacts, Rm, is 
indicated in the labels. 



all contacts with Rn > lOil the rather large mean free 
path of these samples {£ = 80 nm) ensures the fulfillment 
of the conditions for ballistic conduction (see Sect llll l|) 
even if a single contact is hypothesized. Clearly, if several 
parallel contacts are present, they must be necessarily 
ballistic-^. The spectra are normalized, i.e. divided by 
the differential conductance at (being the critical 
temperature of the junction). The experimental curves 
in FiglT3]clearly show the predicted anisotropy^S,, but the 
non-perfect directionality of PCAR prevents the weight 
of the TT-band conductance from assuming the theoretical 
extremal values {w^r = 0.66 for a6-plane tunneling, — 
0.99 for c-axis tunneling). This is particularly clear in 
c-axis contacts, where the single-gap BTK fit (dashed 
line) does not work well and a two-band fit (solid line) 
is instead necessary (with w-^ < 0.99). The values of 
the fitting parameters are indicated in the labels. The 
temperature dependence of the gaps obtained in different 
contacts on single crystals^ is shown in FigfTWb) (left 
triangles). 



2. PCAR in magnetic field 

As mentioned above, the first PCAR measurements in 
MgB2 in the presence of a magnetic field were carried 
out by Szabo et al.^^- in polycrystals. Fig. [H] shows 
the magnetic field dependence of the low-temperature 
PCAR spectra for contacts with a large a6-plane con- 
tribution (a) and a dominant c-axis contribution (b). 
In the first case, the peaks related to are fast de- 
pressed by weak fields and become barely detectable at 
_B = IT, at which the large-gap maxima are still clearly 
visible. In the second case, where no clear peaks related 
to Act are observed in zero field^-, the suppression of 
the TT-gap at 1-1.5 T causes an apparent outward shift of 
the conductance peaks (from about 3 meV to 5 meV in 
Fig fn)) that then start to shrink, because of the suppres- 
sion of the (T-band gap. Actually, the use of polycrystals 
made it impossible to control the direction of both the 
probe current and the magnetic field. This is not ir- 
relevant because of the anisotropy of the critical fields 
in MgB2^i^!^. Indeed, PCAR measurements in single 
crystals^^iS showed that: i) A field of about 1 T "com- 
pletely" suppresses the small gap irrespective of the field 
direction. This does not mean that the tt band becomes 
nonsuperconducting, but simply that above 1 T its con- 
tribution to the Andreev signal becomes experimentally 
undetectable and the conductance curves can be fitted to 
a function like G = Wt^I + {1 — WTr)Ga (where Gjr = 1)- 
Incidentally, this also confirms that the tt band is rather 
isotropic; ii) the direction of the field instead affects the 
behavior of A^-, which is reasonable due to the almost- 
2D character of this band. When B || ab, the A^ peaks 
in the conductance remain clearly distinguishable up to 
9 T, with only some signs of gap closing. Instead, when 
B II c, they merge together at B > 4 T giving rise to 
a broad maximum^^; iii) In any case, at least at 4.2 K, 
Act is very little affected by a field of 1 T, either parallel 
or perpendicular to the ab plane^^; iv) in c-axis contacts, 
the suppression of the 7r-band contribution to the conduc- 
tance at about IT is accompanied by an outward shift of 
the conductance peaks and by an abrupt decrease in the 
amplitude of the spectrum. 

A quantitative study of the effect of the field on the 
gaps requires a fit of the experimental curves. Here the 
main problem is: can the BTK model or its generalized 
version be used to fit the conductance curves when a mag- 
netic field is present? In conventional superconductors, 
Naidyuk et al^ showed that the pair-breaking effect of 
the field can be mimicked, within a generalized BTK 
model, by the broadening parameter P. In other words, 
the total broadening parameter P can be considered as 
the sum of an intrinsic (field-independent) P^ (due to self- 
energy and inelastic scattering effects, see Sect IV 4| and 
an extrinsic P/(-B) due to the magnetic field. This ap- 
proach assumes that the pair-breaking effect of the field 
can be completely represented by the broadening P/(-B) 
while its effect on the DOS is negligible in a first-order 
approximation. 
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Unfortunately, it is not really so. At T = 0, the 
theoretical effects on the DOS curves of both the mag- 
netic field and the actual shape of the FS are very 
relevant^2i22ii2i so that the BTK model (which in its sim- 
plest formulation assumes a spherical FS) even if modi- 
fied to account for lifetime effects and even if a high Z is 
used (tunnel regime), fails in fitting them. One could 
expect something similar to occur in the AR regime. 
However, i) the inadequacy of the BTK model is dra- 
matic only at T = 0, while at low but finite tempera- 
tures (T=4 K) the thermally broadened theoretical DOS 
curves become much more similar to those given by the 
standard BTK-lifetime model; ii) as we will show in the 
following, the model seems to fit much better the experi- 
mental curves than the calculated ones, giving results in 
good agreement with theory. As a matter of fact, the 
attempts to use the BTK model to obtain the field de- 
pendence of the gaps from the conductance curves not 
only have given good results, but have been able to ex- 
tract quantitative information about the diffusivities in 
the two band a^^'^^1^^ . This information is crucial since, 
as shown in Refsj^i^i^iiS, varying the ratio of electron 
diffusivities in the two bands will change the resulting 
macroscopic superconducting properties. Incidentally, 
the possibility to obtain information also on the inter- 
band coupling by means of PCAR measurements in the 
presence of a supercurrent parallel to the interface has 
been recently proposed^^. 

Figure [12] shows the PCAR spectra of afo-plane con- 
tacts on MgB2 single crystals, in magnetic fields parallel 
to the c axis. The curves of panel (a) were obtained in a 
"soft" point contac1j2l, while the curves of panel (b) were 
measured by Naidyuk et al. with a Cu tip^. In both 
panels, the experimental curves (thick lines) are com- 
pared to the two-band generalized BTK fit (thin Hues). 
In panel (a), the fit was carried out with the two-band 
BTK model up to about 1 T, while above this field the 



TT-band features became undetectable so that we took 
cr^ = 1, in agreement with our previous findings in single 
crystal a^^'^'' and with those by Szabo et al. In panel (b), 
instead, the fit was performed with both the a and the tt 
contributions up to the highest field. A similar result was 
also obtained by Bugoslavsky et al. in epitaxial thin films 
of MgB2 , where the 7r-band gap was found to survive up 
to 5 T21. 
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FIG. 15: (a,b) Two examples of magnetic-field dependence 
of afc-plane point-contact spectra in MgB2 single crystals. 
Thick lines are experimental spectra, thin lines in (a) are the 
relevant two-band BTK fit. In (a) the "soft" PCAR tech- 
nique was used, with a Ag-paste contact (from Ref.— ). In 
(b) the contact was made by pressing a Cu tip against the 
sample edge (from Ref.^''). (c) Magnetic field dependence 
of the gaps extracted from the fit of the conductance curves 
in (a) (solid circles) and (b) (open circles). The values of the 
gaps are compared to theoretical predictions in the clean limit 
(solid lines) and in the dirty limitj^ in the particular case 
-Dct = 0.2DTr. The scale for the dashed lines is on the top and 
left axes; the curves for the clean limit have been re-scaled 
vertically and horizontally using the actual Tc of the samples 
and the relevant critical field. 



FIG. 14: Magnetic field dependence of low-temperature nor- 
malized conductance curves measured in MgB2 polycrystals. 
The curves refer to contacts with large aft-plane contribution 
(a) and dominant c-axis contribution (b). Data taken from 
Refs..^ (a) andSi (b). 



The values of the gaps extracted from the fit are shown 
in FigfTsFc) as solid and open symbols. The data are 
compared to the predictions by Dahm et al^ in the 
clean limit (solid lines) as well as to those by Koshelev 
and Golubov^^ in the dirty limit (dashed lines) in the 
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particular case where the difFusivities in the two bands 
are D„ = 0.2£)^, which is the case that allows best fit- 
ting the data in panel (a). The fit of the experimental 
PCAR curves seems to work very well and the resulting 
values of are perfectly compatible with the expected 
field dependence. In all cases^SiZl the broadening param- 
eters of the BTK model, and increase linearly with 
field, giving further support to the distinction between in- 
trinsic lifetime broadening and field-induced broadening 
(proportional to B). 

According to FigfTSTc). a fast suppression of the tt- 
band gap features in weak fields^'^i^i^'^^ indicates that 
the TT-band diffusivity of the samples under study is a 
few times greater than the cr-band one. Since the the- 
oretical field dependence of the gaps in clean limit™ is 
identical to that predicted by the dirty-limit model in 
the case D„ = D-^ (equal diffusivities)^^, the difference 
between our results (a) and those by Naidyuk (b) or 
Bugoslavsky-^ are simply due to sample-dependent vari- 
ations in the diffusivity ratio Du/Dt^, that occurs also in 
different crystals from the same batchii. 

A check of internal consistency of the results described 
above was achieved by studying the partial averaged zero- 
bias density of states (ZBD) in the two bands, N^{Q) 
and iV^(O), whose magnetic-field dependence is predicted 
to depend again on the diffusivity ratic^^. In Refill we 
used the gaps A^- and A^, the weight and the field- 
induced broadening parameters F^''^ obtained from the 
fit of the PCAR spectra in Fig lTST a) to calculate the 
zero-temperature ideal tunneling conductance by setting 
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0. The values of the 



total ZBD as a function of the magnetic field are shown 
in FigfTCTa) (symbols), and compared to the theoretical 
total ZBD N = w^N-^{<S) + (1 - w^)iV<,(0) were 7V^(0) 
and iVcr(O) are those calculated in^ suitably scaled to 
the actual critical field B^2\\c = 6.25 T. It is clear that 
the ZBD follows the theoretical predictions for the case 
Du = 0.2_D7r, in perfect agreement with the conclusions 
drawn from the field dependence of the gaps. 

Bugoslavsky et al^ were able to directly extract the 
partial ZBD from the fit of their PCAR spectra in epi- 
taxial MgB2 films with a suitable model developed and 
tested in conventional type-II superconductors^*'. The 
model is based on the existence of parallel contacts and 
on the fact that a fraction of them (increasing with 
field) occurs in normal-state regions of the sample (vortex 
cores )Z2. Hence, what PCAR measures is an effective av- 
erage over the vortex lattice and the conductance of the 
point contact should be considered as being the sum of 
a "normal channel" and of a "superconducting channel" 
contributions. In MgB2, the normalized conductance in 
magnetic field thus becomes 

G{V) = w^[n^ + {1 ~ n^)G^] + 

+ (l-«;,)K + (l-n,)G,] (32) 

where n^^ and represent the fraction of normal-state 
core excitation and were thus identified^^ with the partial 
ZBD No- and N.^. To reduce the number of free fitting 



parameters, the authors assumed Za- = and used the 
convolution with a Gaussian of width lo to account for 
both the thermal smearing and inelastic interface scat- 
tering. Z and were fixed to their zero-field values. 
The values of and ri^r obtained from the fit of differ- 
ent series of conductance curves, in magnetic fields either 
parallel or perpendicular to the c axis, are reported in 
Fig. [TBf b). The data in i? || c approximately agree with 
the theoretical curves (solid lines) for D„ — 0.5D^. This 
is consistent with the field dependence of the gaps shown 
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FIG. 16: (a,b) Magnetic field dependence of the zero-bias den- 
sity of states evaluated from PCAR. (a) The total ZBD ob- 
tained by simulating the zero-temperature tunneling conduc- 
tance curves (from Ref.^'^) (symbols) compared to the theoret- 
ical prediction in the dirty limit for Da = 0.2_D.^. The perfect 
agreement confirms the result given by the field-dependence 
of the gaps in fig llSf c'). (b) The partial ZBD n,r and ob- 
tained from the fit of the conductance curves with the model 
by Bugoslavsky et al7^ , compared to the theoretical predic- 
tions in the case Da = O.SDtt (sohd lines). The dashed line 
is an estimate of the expected behavior of Na for the case 
H II ab. The results agree with the gap measurements in the 
same samples (shown in (c)) with the persistence of the small 
gap up to 5 T. (d) Magnetic field dependence of the excess 
current lexc (from Ref.— ). Symbols: experimental values of 
lexc from integration of the conductance curves in Fig llSf b') . 
Lines: Iexc{B) calculated from eq. |33]based on the theoretical 
predictions of ref.—. 



21 



in (c), that shows the persistence of the smaU gap up to 
5 T ""^ . Note that the field dependence of A^^ and n.^. does 
not depend on the field direction (indicating that the tt 
band diffusivity is isotropic) while the u-band quantities 
show a marked anisotropy. 

Another interesting way to determine the diffusivity 
ratio D„/Dtj from a PCAR experiment is described in 
Ref73.j where the excess current lexc ~ obtained by inte- 
gration of the reduced dl/dV after subtraction of the 
background - is directly plotted as a function of the 
magnetic field. The result is reported in FigfTCTd) for 
the curves in FigfTSTb). Let us just recall here for 
convenience that the excess current is approximately 
lexc — A/ei?jv; when a magnetic field is applied, a frac- 
tion A^(0, B) of the contact becomes nonsuperconducting 
(vortex cores) and does not contribute any longer to lexc- 
Taking into account the presence of two bands (whose 
partial ZBD behave differently in field) Naidyuk et al. 
arrived to the expression 

Iexc{B) cx w^A^{B)[l- N,{0,B)] + 

+ il-w^)A,{B)[l~N,iO,B)] (33) 

that was used to fit the experimental data. In this func- 
tion, the gaps vales A.^{B) and A(j(i?) as well as the 
ZBD 7V^(0, B) and 7V„(0, B) are taken from the theoret- 
ical curves of Ref<ii suitably scaled to the actual critical 
field. The zero-field values of the gaps and the weight 
were obtained from the fit of the zero-field Andreev spec- 
tra. Fig |16r d) shows that the experimental values of I^xc 
are in good agreement with the theoretical predictions 
in the case — D^^, which further confirms the con- 
clusion drawn from the magnetic-field dependence of the 
gaps in Fig |15f c). It is worth mentioning that both in 
RefjS and2^, different contacts resulted in different val- 
ues of the diffusivity ratio and the curves shown here are 
just an example. 

The low-temperature magnetic-field dependence of the 
gaps in single crystals shown in FigfTSTc) indicates that, 
in these samples, a field of 1 T makes the PCAR spec- 
tra look as if the small gap was completely suppressed, 
but does not seriously affect the cr-band gap. A fit of 
the spectra with a function like G{B = IT) = WttI + 
(1 — WTr)Ga{B = IT) is thus possible and gives a gap 
Act which coincides with the zero-field one - but has a 
smaller uncertainty because of the reduced number of fit- 
ting parameters. If now one subtracts the experimental 
normalized spectrum Gexp{B = IT) from the zero- field 
one, Gexp{B = 0), one obtains a curve that only contains 
the zero-field 7r-band contribution and can thus be fit- 
ted by a 3-parameters function like G{B = 0) — G{B — 
IT) ~ WT;[GTr{B = 0) — 1] from which the small gap A^r 
can be obtained with a small uncertainty. In this fitting 
process, all the parameters must be adjusted to ensure 
a consistency between the different fits. The final result, 
as shown in Refi^, is a rather strict determination of 
the gaps, that turn out to be A^ = 2.80 ± 0.05 meV 
and Ao- = 7.1 ± 0.1 meV, in excellent agreement with 
theoretical predictions of Ref/^. 



The process works well at 4.2 K, and can be extended 
to higher temperatures with some caution, because of 
temperature-dependent anisotropy of the critical fields 



in MgB: 



63.64.65 



. In particular, the field must be parallel 



to the ab plane^i^; in this case the single-band BTK fit 
of the PCAR spectra in a field of 1 T (see Fig fTTf b)) gives 
values of Ao- that agree very well with those determined 
by the two-band fit of the zero-field curves (FigfTTTa^'l 
apart from a much smaller uncertainty. The fit of the 
difference Gexp{B = 0)-Gexp{B = IT) (FigHTJc)) gives 
very good results for the small gap A^ as well. The 
resulting temperature dependence of the gaps^ is shown 
in FigHltd). 
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FIG. 17: (a) Temperature dependence of the zero-field con- 
ductance curves G{B = 0) of a c-axis contact on a MgB2 
single crystal (symbols) and the relevant two-gap BTK fit. 
(b) Same as in (a), but in a field of 1 T applied parallel to the 
ab plane. The G{B = IT) curves are compared to the relevant 
fit with a single-gap BTK model containing only the cr-band 
contribution to the Andreev conductance, (c) Temperature 
dependence of the difference G{B = 0) — G{B — IT), com- 
pared to the relevant single-gap BTK fit that contains only 
the TT-band Andreev refiection conductance, (d) Tempera- 
ture dependence of the gaps and A,r extracted from the 
single-gap BTK fits of the partial a and 7r-band conductances 
(panels (b) and (c), respectively). Data taken from Refs.— 
and^. 

In Refsj^ we also determined by means of PCAR the 
temperature dependence of the critical field i?c2||cj and 



c2\\ab7 



identified with the field that suppresses super- 



conductivity and restores the normal-state conductance. 
For T > 0.8Tc, the critical fields measured by PCAR 
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fall on the curves given by bulk measurements like ther- 
mal conductivitj*^, torque magnetometry25. and specific 
hea^. At lower temperature, they depart from these 
curves and tend to the (larger) values given by resistivity 
measurements^!^. Since the latter rather determine the 
surface critical field -Bc3^, some surface effects clearly 
play a role in PCAR (especially at low temperature). 
Even if PCS is a surface-sensitive technique, this effect 
cannot simply be a surface nucleation of superconductiv- 
ity at a field Bq (5^2 < Bq < Bcs) because the magne- 
toresistivity of the electrode is always much lower than 
that of the MgB2 crystals. 



3. Determination of the electron-phonon spectral function 
by PCS 

Point contact spectroscopy in the normal and super- 
conducting state of MgB2 was also used to obtain the 
electron-phonon spectral function and elucidate the role 
of the in-plane B stretching mode E2g in determining the 
superconducting properties of this compound. In refi^, 
Yanson et al. investigated the point-contact spectra in 
c-axis oriented films, in the superconducting (T < Tc) 
and in the normal {T > Tg) state. They directly mea- 
sured the differential resistance {dV/dl) as well as the 
second derivative of the I-V curves, i.e. cPV/dl^, which 
is proportional to the electron-phonon spectral function 
Q;ppF(a;) (see Sect llll The authors observed clear 
structures in the superconducting state, with a signa- 
ture of the E2g mode of the same amplitude as other 
phonon peaks. Owing to the preferential current injec- 
tion along the c axis, this is in agreement with the calcu- 
lated a^F^(a;)^. However, these structures were found 
to disappear at T > Tc where only much smaller nonlin- 
earities persisted. This indicates a superconducting ori- 
gin of the structures, i.e. due to the energy-dependence 
of the order parameter ("elastic" or self-energy term) 
rather than to the actual inelastic e-ph scattering. This 
point was also addressed theoretically in Ref.^^ where a 
simple asymptotic formula for the order parameter self- 
energy effects in the superconducting point contact was 
derived. Later, the same authors performed PCS mea- 
surements in single crystals in an inverse needle-anvil 
configuration (superconducting crystal as the needle)^. 
The smaller critical field of crystals with respect to films 
allowed at least a partial suppression of superconductiv- 
ity at 4.2K by means of fields as high as 9 T. In low- 
temperature, zero-field spectra with a major contribu- 
tion from a6-plane current injection, the second deriva- 
tive showed a very broad maximum around 60-70 meV 
that was almost insensitive to magnetic field and was 
identified with the signature of the i?2g phonon mode, 
largely smeared by the e-ph coupling, as observed by X- 
ray inelastic scattering^. The d^V/dl"^ spectra of con- 
tacts with a predominant c-axis tunneling contribution 
showed instead much weaker structures such as shallow 
maxima at 30 and 50 mV that were claimed to reflect 



bulk (isotropic) phonons, and were put in connection 
with the first two maxima in the phonon density of states 
or a'^F{u})r^. Further details on this subject can be found 
in refsi^ and^l. 



4- Effect of chemical doping on the gaps of MgB2 

Chemical substitutions in MgB2 were tried very 
soon after the discovery of superconductivity in this 
compounds^. The huge experimental work carried out 
in substituted MgB2 has allowed a deeper understand- 
ing of the pure compound, but has also unveiled a sur- 
prisingly rich an complex physics. Even in the simplest 
effective two-band model^^i^ the quantities needed to 
describe MgB2 are manifold: four Eliasberg functions 
a'^Fij{uj) (where i,j = tr, tt)^, four quasiparticle scat- 
tering rates (intraband, Tu, and interband, Py) two 
densities of states (DOSs), and a prefactor /io to the 
Coulomb pseudopotential (which is a 2x2 matrix whose 
elements only contain the densities of states iVo- and iV^). 
Not all these parameters are independent (for example 
r^CT = (A^^/A^^)r^^, being N^^t^ the zero-bias DOS in the 
two bands) and some of them (like the DOS, the phonon 
spectrum) can be either calculated from first principles 
or determined experimentally. As we will see in the fol- 
lowing, chemical doping in MgB2 always gives a decrease 
in Tc (see Fig fTB)) and in A^,, that can be due either to 
a variation in the DOS (and in the phonon frequencies, 
which however play a minor role) or to an increase in 
interband scattering. Fortunately, the effects of the lat- 
ter are rather easily distinguishable from those of other 
sources of Tc reduction, since an increase in Po-^ sup- 
presses Tc and A^- but increases A^. As a consequence, 
some indications about the effects of doping on the DOS 
and on the interband scattering can be extracted from 
the analysis of the doping dependence of the gaps, mea- 
sured by PCAR, within the two-band Eliashberg theory. 
Indications about the relative role of intraband scatter- 
ing rates in the two bands can instead be obtained from 
the magnetic-field dependence of the ZBD or the ex- 
cess current, as shown in Sect lVII B 21 A detailed review 
about PCAR measurements in doped MgB2 is reported 
in Refi^; in the following we will thus give only a general 
discussion of the main findings. 

In general, a problem with PCAR in doped MgB2 is 
that the structures related to the large gap A^- become 
less and less clear on increasing the dopant content. The 
comparison of the two-band fit to the single-band one 
can clarify whether two gaps are still present or not, but 
sometimes the magnetic-field dependence of the conduc- 
tance curves can be more conclusive in this sense; an out- 
ward shift of the conductance peaks is a strong indication 
in favor of two gaps. Another problem is the determina- 
tion of the actual doping content, especially in the case 
of light atomic species (Li, C) , which makes the trend of 
the gaps vs. the doping content rather uncertain. More- 
over, the doping content is intrinsically inhomogeneous 
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FIG. 18: Bulk critical temperature vs. doping content in vari- 
ous series of substituted MgB2. The doping content for substi- 
tutions in the Mg site (general formula Mgi_a;Ma:B2) must be 
read on the top axis. The content of carbon in Mg(Bi_a;Ci:)2 
is on the bottom axis. The two scales differ by a factor 2 to 
show that all the curves for non-magnetic doping look very 
similar if plotted as a function of the number of substitutional 
atoms per formula unit. The data are taken from Refs.^ (C, 
Ribeiro),-2i (Al, Klein) ,^ (Al, Zambano) (Al, Putti), (Al, 
Mn, Fe, Li, C, Karpinski). 
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FIG. 19; Doping dependence of PCAR spectra in C-doped 
MgB2 polycrystals (a) and single crystals (b). The point con- 
tacts were made with a metallic tip in (a) and with a Ag-paste 
spot in (b). Symbols represent experimental curves, solid lines 
the relevant two-band BTK fit, apart from the top curve in 
(b) that represents instead the single-band BTK fit. Data are 
taken from Ref.-^ (a) and^ (b). 



on the scale probed by PCAR so that different contacts 
on the same sample can even show different gaps and dif- 
ferent . For these reasons, we always prefer to report 
the values of the gaps as a function of the local critical 
temperature . This representation is also the most 
suited to compare the results for different substitutions. 
For example, one can learn that there is an interesting 
universal scaling law of A(j with T^, at least for > 20 
K, independently of the main mechanism of Tc reduction. 

Carbon is the only chemical substitution in the site of 
boron that the structure of MgB2 accepts. PCAR exper- 
iments in Mg(Bi_xC2:)2 were carried out in nearly single- 
phase polycrystals with 0.09 < a; < 0.1 3^°i^^ and in single 
crystals grown at high pressure and high temperature^^ 
with 0.047 < a; < 0.132. The critical temperature de- 
creases on C doping, as shown in Fig[THl although the 
difficulty in determining the actual C content gives rise 
to some minor differences in the actual Tc vs. x curve for 
crystals (black circles) and polycrystals (green circles). 

Fig. [in] shows some examples of normalized PCAR 
curves and the relevant fit at different C contents, in poly- 
crystals and wiresSl and in single crystals^^. Already at 
X = 0.1, the experimental curves seem to be fittable by a 
single-band BTK model with lifetime broadening. Actu- 
ally, at this C content the second gap is still retained^^ 
as it can be easily shown by applying a magnetic field: 
as in pure MgB2, the small gap is fast suppressed and 
the large-gap features emerge clearly^2i2^. As shown in 
Fig fTOT b). in single crystals with the highest doping con- 
tent (x = 0.132) the application of the magnetic field 
makes the conductance peaks decrease in amplitude but 
no shift in energy is observed as long as the applied field 



is much smaller than the critical field (this is certainly 
true for H || ak^^ because C doping increases the critical 
field). The absence of a shift indicates that, if two gaps 
are present, they have very similar amplitude - indeed, 
the two-band BTK fit of all the conductance curves at 
X = 0.132 requires two gap values that are experimen- 
tally indistinguishable. The single-band BTK fit shown 
in FigHHb) gives a gap A = 2.8 ± 0.2 meV. One could 
thus conclude that the "gap merging" is obtained in heav- 
ily C-doped MgB2, also relying on the fact that: i) the 
ratio 2A/kBTc — 3.8 is close to the BCS value; ii) the 
conductance curves recorded in this sample at different 
temperatures all admit a very good single-band BTK fit; 
iii) the temperature dependence of the gap A extracted 
from the fit is perfectly BCS (within the experimental 
uncertainty) and the critical temperature of the junc- 
tion, = 19 K is in perfect agreement with the bulk Tc 
measured by DC zero-field-cooling magnetization. 

However, the situation is not so simple. Fig. [20^ a) 
reports the values of the gaps obtained in C-doped 
crystals^^ and polycrystals^ZiiSSiiifli as a function of the 
critical temperature. It is clear that the two data sets 
agree rather well as far as A^ is concerned, but disagree 
on the values and trend of A^^. In particular, in wires 
and bulk polycrystalline samples there is no tendency to 
the gap merging observed instead in single crystals. To 
analyze the data within the two-band Eliashberg model 
one can use the a- and 7r-band DOS and the phonon 
frequencies calculated for C-doped MgB2^ thus leav- 
ing fiQ and Pcttt as the only adjustable parameters. The 
overall trend of the gaps in C-doped polycrystals is re- 
produced by keeping Fcttt = as in pure MgB2, adjust- 
ing fiQ to reproduce the experimental Tc, and calculat- 
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FIG. 20: (a) Energy gaps in Mg(Bi_a;Ci,)2 as a function of 
the critical temperature, measured in single crystals^* (cir- 
cles) and polycrystals and wires^iiSSiiSl (triangles). Each 
point of the single-crystal series is actually the average of dif- 
ferent gap values measured in different contacts, whose spread 
is indicated by the error bar. (b) Energy gaps measured by 
PCAR in Mgi_a;Ali,B2 as a function of the critical temper- 
ature. Circles are taken fromiS^, squares from^, down tri- 



angles fromiSi and up triangles froro^. In both (a) and (b), 
dash-dot lines represent the gap values calculated within the 
two-band Eliashberg theory by using the DOS and the phonon 
frequencies from ab-initio calculations. Solid lines indicate 
the fit of the gaps vs. Tc obtained by adding to the model an 
adjustable interband scattering rate Tan- 



ing the corresponding gaps (dashed hnes in Fig EOT a)). 
This means that the decrease in Tc, A„ and A^r in these 
samples can be completely explained by band fillingiSl. 
In single crystals, instead, the trend of A.^. can only be 
reproduced by increasing the interband scattering rate 
(solid lines in Fig|20ya)). This contrasts with the theo- 
retical predictioniS^ that substitutions in the B plane (for 
example by carbon) preserving the different parity of a 
and TT bands have little or no effect on Fo-tt- The key to 
this puzzle could be the presence of microscopic defects 
in C-doped single crystals, acting as scattering centers, 
suggested by the doing-induced increase in flux pinning 
and in the normalized resistance^. The nature of these 
defects and the reason why they should be able to create 
interband scattering is however not completely clear—. 
It is worth mentioning that the analysis of the zero-bias 
DOS as a function of the magnetic field in C-doped poly- 



crystals, carried out in refjiSi by using the fitting model 
developed by Bugoslavsky2^, clearly proves that the ra- 
tio D^/Dtt increases from 0.2 (pure MgB2) towards 1 on 
increasing the doping content, indicating that C doping 
(surprisingly) increases the intraband tt scattering more 
than the u-band one. 

Doping in the Mg site has been obtained with different 
chemical species: Al, Li, Mn, Fe. The first two are het- 
erovalent and result in electron and hole doping, respec- 
tively. According to theoretical predictions, Al should 
give the maximum increase in interband scattering (for 
2% of Al a value of Fq-tt = 1.1 meV is predicted, which 
already has measurable effects on the critical tempera- 
ture and on the gapsiS^). On the contrary, Li should 
have little or no effect on To-tt^- PCAR measurements 
in Mgi_2;Ala;B2 polycrystals2^ii2i and crystalsSi up to 
X — 0.2, carried out with either the conventional or the 
"soft" technique, showed an almost linear decrease of A^ 
and A^ as a function of the Al content, in agreement 
with the findings of specific-heat measurements^^. No 
clear tendency of the gaps to merge was observed; an ex- 
tension to higher doping (x = 0.32) was later obtained by 
us in single crystalsiS^. All these results are reported, as 
a function of the critical temperature, in Figl^DTb). As 
in the case of carbon doping, all the data sets agree on 
the behavior of A^, which is also directly related to the 
suppression of the critical field by Al doping22i. Interest- 
ingly, this relationship implies the validity of a clean-limit 
description of the system, as it follows from the analysis 
of the critical field^^iiSi but also from the analysis of the 
zero-bias DOS as a function of the magnetic field (car- 
ried out in ref.^'^'^ by using Bugoslavsky's fitting model*"^) 
which shows that the diffusivity ratio of pure MgB2 is 
preserved on Al doping at least up to x = 0.20. The per- 
sistence of two gaps even at the highest doping content 
is not clear in the spectra; the single-band and two-band 
BTK fit are also very similar to each other although a 
statistical test (the Fisher F test) clearly indicates that 
the latter is preferable for any level of confidence. In any 
case, the magnetic-field dependence of the conductance 
curves shows the outward shift of the conductance max- 
ima on increasing the fiel d^°^'^°^ that always indicates 
the presence of two gaps of different amplitude - even 
though, as explained in Ref.^°^, the suppression of Hc2 
by Al doping prevents the separation of the partial a and 
TT band contributions to the conductance, as we instead 
did in pure MgB2. 

The dependence of A^- on the critical temperature 
shown in Figl^DTbl can be superimposed to that observed 
in C-doped MgEz^^^'^o^ shown in Figl2D)^a). The trend 
of the small gap A^r is similar, but not identical, to that 
observed in C-doped MgB2 polycrystals (see Fig l^Uf a)). 
Actually, a small tendency to an increase in A^ at low 
doping content (with a maximum around = 30 K) is 
very clear, outside the experimental uncertainty, in our 
data on single crystals (sohd circles). Cooley et aZ.^— 
noticed the same trend only in samples produced via a 
long reaction at high temperature so as to reduce the 
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strain and the inhomogeneity in the Al content. This 
would indicate that the enhancement in A^r is intrinsic 
to Al doping but is often masked by other effects, and 
could explain why it is barely detectable in PCAR re- 
sults by the Slovak group^-^iSi (triangles) as well as in 
the results of specific-heat measurements in some poly- 
crystals samples^^iiSi. In any case, this trend cannot be 
reproduced within the two-band Eliashberg model if only 
the proper variation in the DOS^S^ and in the phonon 
frequencies due to Al dopingi^ are taken into account 
(dashed hues in FigEOfb)). Although the effect on the 
DOS is certainly dominantiSl, an increase in interband 
scattering at low Al contents in quantitative agreement 
with expect at ionsiS^ is also necessary to catch the exper- 
imental trend of Ajr in our single crystals. However, in 
order to fit the data, Fo-tt must again decrease for x > 0.1. 
The reason for this is not completely clear but might be 
related to extrinsic effects like inhomogeneity and lattice 
stress (not taken into account in the Eliashberg model) 
that start playing a role for x > 0.1, as suggested in 
j^g£^ji06 g^j-^^ j-,y |-]^g increase in the width of the super- 
conducting transition for x > 0.1-^. 

Mn doping is peculiar for two reasons: i) Mn is ho- 
movalcnt with Mg; ii) its magnetic moment gives rise 
to spin-flip, pair-breaking scattering that is considerably 
larger than the non-magnetic one: indeed, Tc is very fast 
suppressed by small Mn contents (see Fig fTS]) . One can 
thus expect the DOS to be unaffected and the interband 
non-magnetic scattering to play little role in this com- 
pound. PCAR measurements with the soft technique 
were carried out in single crystalsii^ of Mgi_j;Mn2;B2 
with X up to O.OlSii^. The trend of the gaps as a func- 
tion the Tf is shown in FigHHa). For x > 0.004 (i.e. 
for < 33 K) the persistence of two gaps was not evi- 
dent in the conductance curves and had to be proved by 
using a magnetic field. Unfortunately, Mn doping also 
suppresses Hc2 so that, when Tc < 17 K, this procedure 
becomes unreliable and the conclusion that two gaps per- 
sist down to the lowest can be based on: i) the better 
quality of the two-band fitii^, and ii) the fact that the 
presence of a single gap in the low-T^ region would imply 
a sudden change in the slope of the A„ and vs. 
curves^^*^ that is not justified by any observed discontinu- 
ity in the physical properties of the compouncU^. The 
trend of A^ and A^ vs in Mn-doped MgB2 is sur- 
prisingly similar to that observed in Al-doped samples, 
apart from the low-doping enhancement of A^. By the 
way. Act follows the universal scaling law with Tc- Unlike 
in previous cases, the analysis of the data within the two- 
band Eliashberg model can here give precise information 
on the magnetic scattering rates, either interband (F*^ ) 
or intraband (F^). The gap trend can indeed be repro- 
duced very well by using the same phonon spectra, DOS 
values, and Coulomb pseudopotential as in pure MgB2, 
neglecting all non-magnetic scattering rates, and taking 
F^^, F^^ and F*^ as the only adjustable parameters^3£. 
The fit of the gaps vs. Tc indicates a dominant intra- 
band spin-fiip scattering in the a band, F*^, with pos- 
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FIG. 21: (a) Energy gaps measured by PCAR in 
Mgi_iMna;B2 single crystalsii^ as a function of T^. Lines 
indicate the gap trend calculated within the Eliashberg the- 
ory by assuming that Mn doping only gives rise to an increase 
in the a intraband magnetic scattering (r™^.), with a smaller 
contribution from the n — n channel (T^^). (b) Energy gaps 
measured by PCAR (circles^) and specific heat (squares^^i) 
in neutron- irradiated Mg^^B2 polycrystals. Lines indicate the 
trend of the gaps calculated within the two-band Eliashberg 
theory by assuming that: i) the only effect of irradiation is a 
decrease in the cr-band DOS at the Fermi level Na{0) (dash 
lines), and ii) this effect is accompanied by an increase in 
interband scattering Fa-w (solid lines). The straight line rep- 
resents the BCS gap vs. Tc curve. 

sible smaller contributions from either the 7r-intraband 
F*^ or the interband F*^ channels. A large a — a scat- 
tering was predicted theoretically as being due to the hy- 
bridization of the a bands of MgB2 with the d orbitals of 
MitM^,. The dominance of this term on the tt — n or a — n 
channels was instead demonstrated by first-principle cal- 
culations of the electronic structure of MgB2 near a Mn 
impurityi^. 



5. Effects of irradiation on the gaps of MgB2 

The effects of intentional introduction of disorder in 
MgB2 by means of neutron irradiation have been re- 
cently discussed in a review^^^. Here we will just quickly 
mention the results of PCAR measurements in neutron- 
irradiated Mg^^B2 polycrystalsiS,. As explained ir^ii^. 
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the use of isotopically enriched "'^^B, was necessary to en- 
sure a homogeneous distribution of defects in the bulk 
and avoid self-shielding effects. Neutron flux densities 
up to 1.6 X 10"'^^ cm~^s~^) were used, which suppressed 
the bulk Tc down to 8.7 K. The defect distribution is 
very homogeneous, as shown by X-ray diffraction and by 
the small width (0.9 K at most) of the superconducting 
transitionii^. The PCAR measurements were performed 
with the "soft" techniquei^. The severe shortening of 
the electronic mean free patbi^^ made fulfilling the con- 
ditions for ballistic conduction be more and more diffl- 
cult. In the most irradiated sample, even the contact 
with the highest normal-state resistance (40 fi) turned 
out to be in the diffusive regime and showed the typi- 
cal dips a,t V > Vpeakr^, as well as a moderate heating, 
which was shown to be negligible as long as the volt- 
age drop at the junction was of the order of Vpeair^- 
The trend of the gaps as a function of is shown in 
Fig |21f b) (circles). In the region of around 18-19 
K, the results of the two-band BTK fit are shown even 
though a single-band fit (with A « A^) is possible as 
well. The l^a{T^) and A^(T/^) curves clearly indicate 
a transition from two-band to single-band superconduc- 
tivity below 20 K, in excellent agreement with the find- 
ings of specific-heat measurementsiii in the same sam- 
ples (squares). The initial small increase in Ajr suggests 
that neutron irradiation increases interband scattering. 
However, this is not the only effect since a decrease in 
the (T-band DOS (indeed observed experimentally-'^ ^^) is 
necessary as well to approximately reproduce the overall 
trend of the gaps within the two-band Eliashberg model 
(solid line in Fig|2TbV Actually, the DOS decrease is 
dominant and can, alone, qualitatively explain the ex- 
perimental data (dashed lines) and the inclusion of in- 
terband scattering only improves the agreement in the 
high-Tc region. The fit in the low- Tc region (below 20 K) 
is poor but cannot be improved since here both gaps are 
smaller than the BCS value and this is forbidden within 
the 2-band Eliashberg theory (although often observed 
in disordered superconductors^^). 



C. Point-contact spectroscopy in novel Fe-based 
superconductors 

At the beginning of 2008, a new class of Fe-based su- 
perconductors with unexpectedly high critical temper- 
atures - with a record Tc (up to now) of 57 K - was 
discovered . These materials are the first real term of 
comparison for cuprates and thus provide a unique op- 
portunity to test the generality of the theories for high-Tc 
superconductivity and to identify more clearly the con- 
ditions for its occurrence. Many compounds of this class 
have been (and are being) discovered and studied; in the 
following we will only refer to the most widely studied 
families of iron-arsenide superconductors. The so called 
"1111" family includes the compounds REFeAsO (RE= 
rare earth) that become superconducting upon doping 



in the O site, with max Tc=55 K as well as the recently 
discovered oxygen-free Cai_a;REa;FeAsF which shows the 
record Tc = 57 K. The "122" family has general formula 
AFe2As2 (A=Ba, Sr, Ca, Eu) and, upon doping in the A 
site, develops Tc up to 38 K. The state of the research on 
Iron-Pnictide up to May 2009 is (partially) summarized 
in ref-ii^. 

Like cuprates, these materials are layered, with al- 
ternating RE-0 and Fe-As layers, the latter apparently 
playing the key role for the occurrence of supercon- 
ductivity. Bandstructure calculationsii^ and ARPES 
measurements^^ showed that the Fermi surface is quasi- 
2D, and is generally made up of two or three hole-like 
sheets around the F point of the first Brillouin zone, and 
two electron-like cylinders at the M point. This imme- 
diately suggests, in analogy with MgB2, the possibility 
of multigap superconductivity and a dependence of the 
tunneling or PCAR spectra on the direction of current in- 
jection. Experimental indications of multigap supercon- 
ductivity came very soon from measurements of the crit- 
ical field and NQR in LaFeAsOFJ^J^OSS^ from direct 
ARPES measurements in Bao.6Ko.4FeA9i2a, from NMR 
in PrFeAsO0.89F0.il— and so onii^. 

The main issue that PCAR spectroscopy has been 
asked to address is the determination of the number, the 
amplitude and the symmetry of the order parameter(s). 
This information is crucial for the development of theo- 
retical models and to test the existing ones. 

The first PCAR measurements, carried out with the 
conventional needle-anvil technique, seemed to support 
a nodal symmetry of the order parameter, because of the 
systematic observation of a zero-bias conductance peak 
(ZBCP). This happened in LaFeAsOo.gFo.i, where Shan 
et ali^ only observed spectra either featureless or with 
a large ZBCP and much smaller additional features. By 
fitting the conductance curves with the 2D BTK modePi 
(see Sects IV 21 and IV 6p . they obtained a c?-wave gap 
A = 3.9 ± 0.7 meV, corresponding to 2A/fc_BTc = 4.11. 
Also Wang et al}^^ often observed a ZBCP in polycrys- 
talline SmFeAsOo.gFo.i and interpreted it as a signature 
of nodal gap. All their spectra admitted a fit with the 
(i-wave 2D BTK model, including the few spectra with 
no ZBCP. 

The interpretation of the ZBCP as an indication of 
nodal symmetry was however soon questioned. Yates et 
al^ performed PCAR in oxygen-deficient NdFeAsOo,85 
with Tc = 45.5 K and observed that a ZBCP (always 
vanishing at T^) develops on increasing the pressure ap- 
plied by the tip (and thus on decreasing Z and Rm)- 
This made the authors warn that the ZBCP may be 
an artifact, not related to the gap symmetry. A tem- 
perature dependence of the ZBCP incompatible with a 
d-wave symmetry was found by the same authors also 
in TbFeAsOo.gFo.i with Tc = 50 K, where the pres- 
ence or absence of the ZBCP was found to even de- 
pend on the sample region probed by PCARf^. Samuely 
et alr^ observed, in NdFeAsOo.gFo.i, a large predomi- 
nance of low-temperature spectra without ZBCP. How- 
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ever, the ZBCP was found to emerge at a temperature 
T* cx to grow in amplitude until it overwhelms 

the gap structures (as ini^l)^ and finally to disappear 
at T/. A ZBCP very robust against the magnetic field 
was found in low-resistance junctions by Chen et al}^^ 
who performed PCAR measurements in polycrystalline 
SmFeAsO0.85F0.15 with bulk Tc = 42 K. In contrast with 
ref^l^^, Chen et al. also found that the ZBCP is very lit- 
tle affected by magnetic fields up to 9T - which excludes 
its relationship with a d-wave symmetry of the gap. 

The actual nature of this peak is not completely clear. 
Kondo scattering coming from the magnetic impurities 
in or near the barrier-i2^, or to the magnetic moment of 
Nd or Tbi21ii2^, can be excluded. The same holds for 
intergrain Josephson couplingi^. Owing to the decrease 
in i?jv that accompanies the increase in tip pressure, one 
could hypothesize that the ZBCP is due to the contact 
not being in the ballistic regime so that critical current 
effects^ or heating in the contact region occurs. It is true 
that, as shown in Ref.-^, the ZBCP coexists with clear 
gap features at finite energyi^ii21ii2^, but this does not 
necessarily exclude the possibility of parallel ballistic and 
thermal contacts. What instead seems to exclude this 
picture is that soft point-contact measurements carried 
out in LaFeAsOi_^F:r^ii23 and SmFeAsOo.gFo.?^'^ 
never gave evidence of ZBCP irrespective of the contact 
resistance, which ranged between a few Q, and more than 
250 Oi^. Therefore the ZBCP is probably related to the 
pressure rather than to the contact resistance, and this 
points towards its relationship with local lattice defor- 
mations. 

Apart from the zero-bias anomaly, the spectra mea- 
sured by Yates et al. in NdFeAsOo.85 show conductance 
maxima in the low-pressure ZBCP-free spectra (bottom 
curve in Figl^H^a)) which evolve smoothly into shoul- 
ders in the high-pressure ones (top curve in Fig lSST a)). 
The position of these features is robust against pressure 
and, if taken as an indicator of the gap amplitude, it 
gives A = 7 meV corresponding to 2A/A;b Tc — 3.57. 
A BTK fit (in s wave-symmetry) of the bottom spec- 
trum in Figl^Ha) gives A ~ 6 meV. In F-doped Nd-1111 
(NdFeAsOo.gFo.i)— , the low-temperature spectra are 
sometimes featureless (bottom curve in Fig|22rb)) but 
more often display conductance peaks at = ±5 — 7 mV 
and shoulders at ^ = 10 mV (top curve in Figl^Wb)). 
which makes them strikingly resemble the PCAR spectra 
in MgB2. Indeed, they were rather well fitted by a two- 
band BTK model where G = aGi + {\ — a)G2 being G the 
normalized junction conductance, and a the (unknown) 
weight of the contribution of band 1. The best-fitting 
values of the gaps are Ai ~ 4 — 6 meV and A2 ~ 9 — 13 
meV that, being Tc = 45 K, imply gap ratios of 2.6 and 
5.7, respectively. Multigap features were clearly observed 
in Tb-lllliiS^, in the form of peaks at about 4.5 and 8 
meV (see top curve in FigllKc)). The two- band BTK fit 
of the spectra gave, at low temperature, Ai = 5 — 6 meV 
and A2 = 8 — 9 meV. The PCAR measurements carried 
out by Wang et al. in polycrystalline SmFeAsOo.gFo.i 



gave two distinct families of spectra, with conductance 
maxima in different positions. Two of them, with no 
ZBCP, are shown in Figl^DJd). Their fit with a single- 
band d-wave model gave Ai = 10.5 ± 0.5 meV and 
A2 — 3.7 ±0.4 meV^^^. Similar gap values were given by 
the two-band d-wave fit of the only spectrum with clear 
two-gap structures (and a ZBCP). However, the spectra 
shown in Figl^Hd) also admit a two-band BTK fit in s- 
wave. Although the very small amplitude of the signal 
requires huge values of the broadening parameters, the 
two s-wave gaps Ai ~ 5 meV and A2 = 14— 18 meV can 
be obtained, which are much higher than those given by 
the single-band d-wave fit. 
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FIG. 22: A few examples of normalized low-temperature 
PCAR spectra measured in various Fe-As compounds of the 
1111 family (indicated in the labels) by means of the con- 
ventional needle-anvil technique. The data are taken from 
Refs.i2I (a),i22 {h)2^ (c)^ {d)^ (e). Some spectra are 
vertically offset for clarity. In (c), the bottom curve is mea- 
sured in a different region of the sample. 
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The high-resistance spectra measured by Chen et al. 
in Sm-1111, showing no ZBCP up to T^, can be fitted in 
the whole temperature range with a s-wave BTK model. 
The raw conductance curves reported in Refsfi2£ii2^ fea- 
ture a characteristic left-right asymmetry ~ also observed 
in Nd-1111 by Samuely et al^ and by us in La-lllLi^ 
and Sm-llllii2i - probably intrinsic to these materials, 
as it was in the case of cuprates^. The normalized spectra 
(of which three examples are reported in Figl22je)) have 
very high Andreev signal, always show clear conductance 
peaks at F ~ ±7 meV related to a superconducting gap, 
and often also display additional structures at higher en- 
ergy that were seen by the authors as extrinsic features 
and thus disregarded in the fitting procedure. The single- 
band BTK fit done in^^° clearly reproduces only the two 
peaks at about 7 nieV and, in the best cases, the zero- 
bias minimum between them. The resulting gap is very 
robust against the contact resistance and turns out to be 
almost perfectly BCS, i.e. A 6.67 ± 0.15 meV at low 
T, corresponding to 2A/kBTc ~ 3.68. 

The existence of a single BCS gap in iron pnictides is 
actually surprising since, at the present state of knowl- 
edge, there seems to be no possible weak-coupling mech- 
anism able to justify the high Tc of these compounds. 
The electron-phonon coupling is very weaki^ while the 
coupling mechanisms mediated by spin fluctuationsi^l 
proposed for these materials require a strong interband 
coupling, and are rather unlikely to give the same BCS 
gap on all the sheets of the Fermi surface. On the other 
hand, many results, not only from PCAR, speak in favor 
of a multiband picture. The same spectra measured by 
Chen et al. may also give indications of multiple gaps, 
if the additional shoulders at > Vpeak are interpreted 
as the hallmarks of a larger superconducting gap. This 
interpretation is questioned by Chen et al. because the 
position and the amplitude of the additional features are 
contact-dependenti^^. This argument holds in conven- 
tional superconductors but, for example, it does not in 
MgB2 (where the observability of the Ao- peaks actually 
depends on the current direction) or in cuprates (where 
a large spatial inhomogeneity in the gap values has been 
observed by STMi^S). 

The differential conductance curves measured in "soft" 
point contacts on SmFeAsQi-j^F^r polycrystals synthe- 
sized at high pressure with x ~ 0.20 {Tc — 52 K) and 
X = 0.09 [Tc — 42 K)i2i generally look very similar to 
those by Chen et al. and always present two clear con- 
ductance peaks and additional shoulders that we chose 
not to disregard in the fit. Similar structures, even more 
marked, were observed by soft PCAR in LaFeAsOi-^^Fa, 
polycrystals with bulk Tc = 27 K. The normalization of 
the spectra prior to fitting is complicated by the fact that: 
i) the normal-state spectrum at low T is not accessible 
due to the high Hc^r^ ii) the normal state at Tc is not flat 
but shows a zero-bias hump in Sm-1111 and a pseudo- 
gaplike feature in La-1111, both progressively washed out 
on increasing temperature until, around the Neel tem- 
perature of the parent compound (about 140-150 K), the 







0.054 








O 








a 


0.052 


o 






■D 




0.050 








O 

o 




0.048 




^ « 
CO tj 

o ^ 
c o 
o 



CD 
O 

B a 

o 

13 ^ 
C ■ — - 
O 
O 



(D 

"O o 
CD c 
^ CC 
CO tj 

o ^ 
c o 
o 



1.15 
1.10 
1.05 
1.00 



(b) p 


LaFeAsO„/„ 


Tj,*=28.6 K /->. / 




A =2.75 meV 


A,=7.9 meV 



-30 -20 -10 
0.08 



0.05 



10 20 30 



0.04 



1.20 
1.15 
1.10 




(d) A 


A SmFeAsO„/„^ ■ 


T "=52 K ^ 

c X/ 




— ij=5.7 meV 


A,=1 9.0 meN/N'-^ars^ 







1.05+ 

-60 -40 -20 20 40 60 



Voltage (mV) 



FIG. 23: (a) The raw conductance curve in a soft point con- 
tact on La-1111 measured at 4.3 K (red curve) and in the 
normal state at Tc (black curve), (b) The normalized low- 
temperature curve (red thick line) and the relevant 2-band 
generalized BTK fit (thin blue line). (c,d) Same as in (a) and 
(b) but for Sm-1111. Data are taken from Refs.^^^ and^^. 



conductance becomes fiat (but retains its right-left asym- 
metry) . This suggests that the normal state conductance 
may change with temperature also below Tc. In the case 
of Sm-1111, we chose to divide all curves at T < Tc 
by the normal state at Tc. In La-1111, we tried differ- 
ent normalizationsi^ and showed that their choice has 
a small effect on the small gap (less than 2%), while it 
can change the larger one by about 10% though preserv- 
ing its trend as a function of temperature, magnetic field 
and critical temperature. Fig. [23] reports an example 
of raw conductance curves at low temperature (4.3 K) 
and at T^ in La-1111 (a) and Sm-1111 (c). The low- 
temperature normalized curves are shown, together with 
the relevant two-band BTK fit, in panels (b) and (d). 
The resulting values of the nodeless gaps are reported in 
the labels. At least in Sm-1111, the normalized curves 
still feature an asymmetry that might be due to the nor- 
malization (i.e. the asymmetry of the normal state might 
depend on temperature) or might be an intrinsic feature 
of these compounds, as it was for cuprates. Certainly, it 
increases the uncertainty on the gaps (A2 in particular) 
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since the model can only fit either side of the experimen- 
tal spectrum (in this case, the right-hand one). 

The temperature dependence of the gaps in La-1111 
and Sm-llll"'^i"3ii34 jg reported in FigllUa) and (b). 
In La-1111 (a), three curves are shown, obtained in con- 
tacts with different and thus possibly made in regions 
with different local doping. Note that = 31K corre- 
sponds to the very beginning of the resistive transition. 
At low temperature, the large gap A2 decreases on in- 
creasing and apparently disappears at = 31 K, 
while the small gap Ai increases. This finding could be 
compatible with recent predictions of a doping depen- 
dence of the gapsi^Sii^ although any conclusion in this 
sense is definitely premature. The anomalous T depen- 
dence of the gaps in La-1111, where A2 seems to close 
at a T* < above which Ai shows a "tail" , could be 
due to the shortening of the mean free path on increasing 
temperature, so that the junction ceases to be ballistic at 
a voltage that decreases with temperature. This would 
mean that the superconducting features are progressively 
weakened starting from the high-bias ones. An indication 
in this sense is given by the high-temperature curves that 
are more peaked at zero bias than expected, but if this 
was the case an apparent decrease in with respect to 
Tc could be expected, which is instead not observed. 

FigEHb) shows the gaps values extracted from the fit 
of various soft-PCAR curves in Sm-lllli^2ii2i. Here the 
situation is clearer: Ai almost coincides with that de- 
termined by Chen et uIf^ (although in our case = 
51 - 53 K while ii>i^ Tc = 42 K); its values are well 
reproducible and follow a BCS-like temperature depen- 
dence up to T^. The trend of A2 vs. T in each set 
of data is compatible with a BCS-like curve, but the 
absolute values are scattered within a region bounded 
by two BCS-like curves with 21^2/kBTc = 7-9. This 
spread could both due to the residual asymmetry of the 
normalized conductance curves and to the uncertainty 
introduced by the normalization. The possibility has 
also been explored theoretically that two large gaps exist, 
of similar amplitude and thus virtually indistinguishable 
by PCARiiSiiil. The gaps obtained by fitting Wang's 
curves in Figl^^c) (green circles) with a two-band s- 
wave 2D BTK model turn out to be in agreement with 
the other results. 

In summary, the interpretation of PCAR spectra in 
1111 compounds within a multigap picture gives reason- 
able results, with a certain degree of universality for the 
different compounds: i) the low-temperature small gap is 
a little smaller than BCS in amplitude: 2Ai/A:bTc ranges 
between 2.23 and 3.44 in La-llll^^^, between 2.54 and 
2.95 in Sm-lllli34 and is around 2.1 in Tb-lllli^S; ii) 
the large gap is larger or much larger than BCS with 
2A2/A;bTc equal to about 4 in Tb-1111, ranging from 
6.42 to 8.68 in La-1111 (when A2 is detectable) and 
from 6.7 to almost 9 in Sm-1111. Such high values are 
confirmed, among others, by ARPES measurements in 
Nd-llllii? and by infrared ellipsometry in Sm-1111^'^'^. 
They indicate that a non-conventional paring mechanism 
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FIG. 24: Upper panel: Energy gaps in La-1111 as a func- 
tion of normalized temperature. The three sets of curves re- 
fer to different contacts with different (indicated in the 
legend)^^. Lower panel: The gaps in Sm-1111 as a func- 
tion of the normalized temperature^^. The data taken from 
Ref.i^ (black circles) are also included, as well as the result 
of the two-band s-wave fit of the low-temperature curves by 
Wang et alM^ shown in Fig. [22jc) (big green circles). 



is taking place - and indeed they can be obtained within 
the Eliashberg theory^**^- by supposing a spin-fluctuation- 
mediated pairing mechanism related to the nesting of the 
Fermi surfacei^ii^ and which should give rise to the so- 
called s± superconductivity, with nodeless order param- 
eters of opposite sign on the electron-like and hole-like 
sheets of the Fermi surface. Unfortunately, PCAR in 
polycrystals cannot give indications about this expected 
TT-phase change. A definite answer to the single gap vs. 
multigap debate could come from ARPES, while phase- 
sensitive techniques are needed to establish whether the 
symmetry is really s±, but probably these developments 
will need the growth of large enough single crystals. It 
must be said that preliminary measurements in a5-plane 
contacts on Sm-1111 single crystals ^"^^ have perfectly con- 
firmed the results mentioned heroic, supporting the ex- 
istence of two nodeless gaps. 

The 122 compounds have been soon grown in the 
form of large crystals, with markedly layered struc- 
ture and easily cleavable. ARPES measurements^^ in 
(Ba,K)Fe2As2 gave unambiguous evidences of multiple 
nodeless gaps (A2 ~ 12 meV on the two small holelike 
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FIG. 25: (a) An example of spectra at low temperature (red 
line) and at (black line) measured in afe-plane contacts in 
hole-doped Ba-122. (b) The normalized curve (red line) and 
he relevant two-gap BTK model (blue line), (c) Temperature 
dependence of the gaps, (d) Typical spectra measured in c 
axis contacts, at low temperature (4.5 K), just above Tc (30 
K) and well above Tc (60 K). Al the data are adapted from 
Refs.iia andi^. 



and electron-like FS sheets, and Ai ~ 6 meV on the large 
hole-like FS). Directional PCAR measurements with a 
Ft tip were performed in hole-doped Bao.55Ko.45Fe2As2 
single crystals by the Slovak group The majority 
of a6-plane point contacts showed a broadened Andreev- 
reflection feature at zero bias (but no ZBCP) and a pseu- 
dogaplike feature in the normal state that persists well 
above and is very similar to that observed in La- 
llllJ^. As in La-1111, it is not clear whether the con- 
temporaneous observation of superconducting signatures 
and a pseudogaplike feature means that they spatially 
coexist or, instead, they belong to spatially separated 
superconducting and antiferromagnetic phasesii^. 

In some contacts, clear double-gap structures were ob- 
served, i.e. symmetric peaks at ±2 — 4 meV and addi- 
tional shoulders at higher energy (about 10 meV) as in 
Fig lSsT a). In these contacts, the normal-state spectrum 
at showed a broad hump at zero bias (black curve in 
Fig l25f a)). The conductance curves were then normal- 



ized to the normal-state spectrum and fitted with success 
to a two-band BTK model, as shown in Fig EST b). The 
low-temperature gaps obtained in several contacts are 
Al = 2.5 — 4 meV and A2 = 9 — 10 meV, corresponding 
to gap ratios ^I^yI^bTc = 2.5-4 and 2Ai/fci3rc = 9-10. 
Both Al and A2 were found to follow rather well a BCS- 
like trend and to close at the same (see Fig l^Sl c)). 
The values of these gaps differ from those observed by 
ARPES in Bao.6Ko.4Fe2As2i20 Ai ~ 6 meV and A2 ~ 12 
meV but this difference is partly justified by the different 
Tc of the samples (T^ = 37 K in Ref J^, = 23 K in 
Ref^). 

No trace of gap features were instead observed in c- 
axis contacts (Fig ESf d)). but only a V-shaped conduc- 
tance progressively flattening on increasing temperature, 
with no apparent signature of the crossing. The filling 
effect cannot be simply explained by the thermal broad- 
ening of the spectra and continues up to about 70 — 80 
K, the temperature at which the magnetic transition in 
the system takes place^^. 

This marked anisotropy of the spectra is interesting 
and it is tempting to associate it with an anisotropy 
of the FS. However, the 122 compounds seem to have 
a nearly-3D FgisSiiis gp the complete absence of 
gap signature along the c direction is difficult to explain; 
moreover, as the authors discuss, the inability to observe 
the gaps in c axis contacts might be as well related to 
surface contamination or reconstruction. 

PCAR measurements were also carried out in the 
electron-doped system Ba(Feo.93Coo.o7)2As2 with bulk 
Tc = 23Ki^. No multigap features were ever observed 
in this material; the spectra allowed instead a fit to 
a standard BTK model, although with large broaden- 
ing, that gave a single isotropic gap of about 5-6 meV. 
However, ARPES measurements in the same system 
Ba(Feo.925Coo.o75)2As2 (Tc = 25.5 K) gave evidence of 
two nodeless gaps Ai = 4.5 meV and A2 = 6.7 meV-i^. 
Evidences of multiple gaps have also been observed in re- 
cent (preliminary) soft PCAR measurements in the same 
material. 



D. PCS and PCAR in borocarbides 

Since their discovery in 1994, borocarbides with for- 
mula i?T2B2C2 (i?=rare earth, T=transition metal, usu- 
ally Ni) have been the subject of an intensive study. 
Their crystal structure resembles that of high- Tc cuprates 
but their electronic properties are rather those of 3D met- 
als. Point contact measurements have been performed 
in these materials both to determine the electron-boson 
spectral function (which allows clarifying the nature of 
the superconducting coupling) and to investigate the 
symmetry of the superconducting state and/or the pos- 
sibility of multiband superconductivity. 

Among the non-magnetic quaternary borocarbides, 
YNi2B2C is probably the most widely studied. Initially 
claimed to be a s-wave superconductor, it has later been 
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found to display a large anisotropy in the superconduct- 
ing state but basically isotropic normal-state properties, 
which suggests an anisotropic order parameter. The pro- 
posed symmetries were s + g (or anisotropic s-wave), 
with point nodes in the ab plane (and precisely along 
the [100] and [010] directions) and the c?-wave symmetry, 
with nodes both in the ab plane (in the [100] and [010] di- 
rection) and along the c axis (where instead the s + g gap 
has a finite value). However, already in 1998, well before 
the discovery of MgB2 , Shulga et alM^ proposed a two- 
band model to explain the temperature dependence of the 
critical fields measured in YNi2B2C and LuNi2B2C. As 
we will see in the following, the debate about the struc- 
ture and the amplitude of the OP in these compounds 
is far from being settled, but the multiband scenario is 
definitely the most supported by the many experimental 
results. 
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FIG. 26: (a) Temperature dependence of the gaps in 
YNi2B2C single crystals along different directions, measured 
by Naidyuk et alM^ (open symbols) and by Raychaudhuri et 
alw^ (solid symbols). The measurements with / || a down to 
300 mK are unpublished (by courtesy of P. Raychaudhuri). 
Solid lines are BCS-like temperature dependencies, (b) Tem- 
perature dependence of the gaps in LuNi2B2C single crystals 
determined from the two-band fit of the PCAR conductance 
curves by Naidyuk et al^. (c) Gaps in ErNi2B2C given 
by the two-band fit of the PCAR spectra along the ab plane 
(squares) and along the c axis (circles) taken from Ref.'^^--. 

Low-T PCS measurements in the normal state of 
YNi2B2C, obtained by suppressing superconductivity 
with a magnetic fieldi^ii^ identified a well-resolved 



maximum in the second derivative cPV/dP at 12 meV 
that corresponds to a soft phonon mode contributing to 
about 90% of the total electron-phonon coupling. Other 
phonon peaks at 20, 24 and 32 meV were not observed. 
This indicates a superconducting coupling mainly me- 
diated by soft phonons, which usually makes an un- 
conventional gap symmetry rather unlikely. Gap mea- 
surements in the superconducting state were carried out 
by PCAR in single crystalsi^^iiSiiS and c-axis oriented 
films'-—. In single crystals with Tc = 14.6 K Raychaud- 
huri et al. performed directional PCAR measurements 
by injecting the current either along the a or the c axis. 
Their conductance curves always admitted a single-band 
BTK fit that however gave clearly different gaps, i.e. 
a small A/||„ = 0.37 — 0.49 meV and a much larger 
A/||c = 1.8 - 2.2 meV. The ratio F/A was fou nd to 
be larger for / ]| a, possibly indicating (see Sect IV 4)) 
a greater angular variation of the gap for a6-plane con- 
tacts, as expected for a s + g symmetryi^. The tem- 
perature dependence of the large gap, which closes at 
= Tc = 14.6 K, was found to slightly deviate from 
a BCS-like curve (solid circles in Fig l^Sf a)). which in 
principle is compatible with a gap with nodes. However, 
the small gap A/||a was found to be fast suppressed on 
heating, falling below a BCS-like curve with Tc — 14.6 
K (dashed line in Fig ESf a)) and to become undetectable 
above 8 K. The recent measurements of A/||a down to 300 
mK carried out by the same group and shown in Fig |26r a) 
(solid squares) indicate that, for T < 3.5 K, this gap fol- 
lows very well a BCS-like curve with Tc = 4.75. As re- 
marked by the authorsiS, this temperature dependence 
is inconsistent with a, s + g symmetry or with any gap 
function of the form A{k) = Ao/(fc), and is indeed more 
compatible with a picture in which the two gaps open 
on different, weakly coupled bands^ (see FigfTTj). The 
magnetic-field dependence gives similar results^^, i.e. 
the small gap A/n^ij is fast suppressed by the magnetic 
field and "disappears" well below H^2- This is very simi- 
lar to what happens in MgB2 when the diffusivities in the 
two bands are differentil. Also the zero-bias DOS, cal- 
culated by using N{E) = + iV) / {E + iT)^ - A2] 
and taking E = 0, increases with field^^^ in a way sim- 
ilar to that predicted by the two-band model for dirty 
superconductors^i. If this is the case, the fact that the 
small gap almost only contributes to the conductance for 
/ ]| a suggests that it opens on the nearly-cylindrical, 
fast-electron Fermi surface sheets that have the maxi- 
mum cross-section (and dominate the conductance) for 
I \\ a but play almost no role for I \\ oiS. The multi- 
band picture is also strongly supported by the effect of 
nonmagnetic (Ft) doping on the critical temperature and 
the upper critical ficld^^^. The decrease of both Tc and 
Hc2 (the latter for either orientations of the field, H ]| a 
and H ]| c) and their subsequent saturation on increas- 
ing doping was indeed shown to be explainable, within 
a two-band picture, as being due to an increase in in- 
terband scattering as in doped MgB^i^. Going back to 
PCAR measurements, further support to the multiband 
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picture also came from the magnetic-field dependence of 
the excess current I^xc obtained by integrating the nor- 
malized PCAR spectra in YNi2B2C c-oriented filmsi^ 
and single crystalsi^ii^. As in MgB2 (see FigfTOTd)) 
Iexc{B) shows a positive (although small) curvature. In 
these measurements, however, the small gap measured 
by Raychaudhuri et al. was never seen. The conduc- 
tance curves were found to display only one peak (at 
either positive or negative bias) and were thus fitted to 
a single-band BTK model giving a distribution of gap 
values. In single crystals with = 15.4 K Naidyuk et 
al}^^ found for the a axis A[ioo] = 1.5—1.7 mcV, for the 
c axis A[ooi] = 1.8— 2.5 meV, and for the [110] direction 
A[iio] = 1-0 — 2.5 mcV. A representative temperature 
dependence of the gaps for / || a and / || 6 is shown in 
Fig ECT a) as open symbols. 

In LuNi2B2C single crystals, PCAR measurements car- 
ried out by Bobrov et alr^ gave spectra with a single 
conductance peak that were however shown to be poorly 
fitted by a single-band BTK model, and a little better, 
but still unsatisfactorily, by a two-band one. Two differ- 
ent continuous distributions of gap values were thus used 
to reproduce the shape of the low-temperature spectra, 
with A/||Qb ranging between 1 meV to 3.35 meV (with 
maxima at 2 meV and 3.1 meV) and A/||c ranging from 
0.7 to 4 meV (with maxima at 1.8-2 meV and 2.5 meV). 
j ^i54,i6i ^ a more conventional two-gap approach was used 
to fit the data in either direction, although the spec- 
tra always featured a single conductance peak - possibly 
because the two gaps are too close to each other to be 
clearly distinguishable. As shown in Fig lSBT b). a qualita- 
tively similar (but quantitatively different) gap trend as a 
function of temperature was observed for I \\ ab (squares) 
and / II c (circles). The results seem to indicate the ex- 
istence of two bands with a weak, anisotropic interband 
coupling^^. In a more recent PCAR study in LuNi2B2C 
single crystals '^^^ the spectra were fitted to a single-band 
BTK model and gave almost equal gaps in the [DDI] (c 
axis, A[ooi] — 2.4 meV) and [110] directions (Ajug] — 2.6 
meV), which clearly excludes a s + g symmetry. Both 
these gaps approximately follow a BCS-like temperature 
dependence with 2A/kBTc equal to 3.4 and 3.6, respec- 
tively - although an upward deviation of the experimen- 
tal points was observed at low temperature, as in films 
Ref These results may be compatible with those by 
Naidyuk et alJ^ in the ab plane, since the single-band 
fit can reasonably give an "average" gap with respect to 
the two-band fit and is probably only weakly sensitive to 
the faster depression of the small gap. Along the a axis, 
a rather wide distribution of gap values Ajj^qq] — 1.6 — 2.5 
meV was observed, all disappearing at the bulk Tc - sur- 
prisingly similar to the findings by Bashlakov et al.^^^ in 
YNi2B2C films. 

Various other compounds of the family were studied 
by PCS in the normal state and PCAR, especially by 
the Ukrainian group, and it would be impossible to ac- 
count here for all their results. Let us just briefly men- 
tion for its interest the magnetic compound IIoNi2B2C 
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FIG. 27: Point-contact electron-boson spectral function 
a%CF{eV) of HoNi2B2C (thick line) compared to the neu- 
tron phonon DOS of the isostructural compound LuNi2B2C 
(thin line). The correspondence of the structures at 3, 10, 16, 
22 and 50 mV is clear. The additional peaks due to magnetic 
and crystal-field excitations are also indicated (from Ref.-iSi.) 

with a low-temperature commensurate AFM state whose 
Neel temperature = 5.3 K is smaller than Tc = 8.5 
K. PCS measurements of a'ji^F{uj) allowed identifying 
structures at 16, 22, 34 and 50 mcV corresponding to 
peaks in the phonon DOS of isostructural LuNi2B2C, 
plus a peak at 3 meV related to the magnetic ordered 
state and indeed disappearing at Tat, and a peak at 10 
meV possibly due to a coupling of carriers to crystal-field 
excitationsi^>i^ (see Figl27]). PCAR spectra in the su- 
perconducting state first of all showed a negligible depen- 
dence on the current direction and admitted single-band 
BTK fit; the low-temperature gap is Aq = 0.95 meV, and 
decreases on increasing temperature following a BCS-like 
curve but disappears at T* — 5.6 K, well below Tc- Be- 
tween T* and Tc, an unconventional gapless state was 
proposed, possibly due to the peculiar spiral magnetic 
order in this temperature range. Below T/v, the super- 
conducting state coexists with an antiferromagnetic or- 
der. To explain this phenomenon, a separation of the two 
phases in the k space has been proposed. In particular, 
superconductivity below Tn should survive only on a sin- 
gle (isotropic) sheet of the FS with no contributions from 
the Ho 5d states. This picture has been directly put in 
connection with the results of critical field and anisotropy 
measurements as a function of temperature in the same 
compoundi^. 

More recently, PCAR measurements have been per- 
formed in the ErNi2B2C compound^^^, with Tc = 11 
K and a low-temperature incommensurate antiferromag- 
netic order with spin density wave between 2 K and 
Tn ~ 6 K. To account for the magnetic pair-breaking 
effect, a suitable model by Beloborod'koi^ was used to 
fit the spectra. The model contains as adjustable pa- 
rameters A (order parameter, OP), Z (barrier height) 
and 7 (magnetic scattering rate). The gap Aq is re- 
lated to the OP by the formula Aq = A(l - 72/3)3/2265^ 
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The single-band fit of the PCAR spectra as a function 
of temperature gave A(T = 2K) ~ 1.8 meV for both 
I \\ ab and I \\ c. A two-band fit was also carried out be- 
cause of the claimed unsatisfactory quality of the single- 
band one. This fit gives two OPs Ai « 2 meV and 
A2 « 1.2 meV whose temperature dependence is shown 
in Figl^STcV The anisotropy in this case is small, but a 
very unconventional behavior is observed because of the 
AFM order below T^. The observation of one or two OPs 
at low temperature clearly proves the coexistence of su- 
perconductivity and AFM order, also observed by laser 
photoemission. On heating, Ai and A2 first increase 
(possibly because of the weakening of the AFM order), 
which is consistent with previous findings by tunnel and 
laser photoemission spectroscopy as well as with the pre- 
diction of some theories of coexistence of superconductiv- 
ity and antiferromagnetic state^^. In the paramagnetic 
state above 6 K, the OPs follow BCS-like curves (whose 
extrapolation would give T^^i ~ 11.3 K and T^'g ~ 14.5 
K) and finally abruptly disappear at Tc = 11 K. The ex- 
trapolated values give for both OPs a ratio 2A/A:bTc 
in the strong-coupling regime. To fit the curves, also 
the weight of the larger OP A2 had to be varied in a 
non-monotonic way, although it generally decreases on 
heating. The magnetic scattering rates also evolve non- 
trivially as a function of temperature, and 72 is always 
greater than 71. This results were interpreted as indi- 
cating that: i) different bands are differently affected by 
magnetic order; ii) the part of the FS that develops the 
larger OP A2 tends to diminish on approaching Tc] iii) 
there is a FS separation with distinct superconducting 
and magnetic bands (or FS sheets); a more detailed anal- 
ysis of the results in Ref.^^^ indicates that approximately 
half of the FS is nonsuperconducting. 

It follows from the above that the result of PCAR 
measurements in borocarbides carried out by different 
groups are often in disagreement with one another. More- 
over, the coexistence of magnetic orders of some kind 
and superconductivity in some of this compounds pre- 
vents any tentative description of this class of com- 
pounds as a whole. Nevertheless, various hints strongly 
suggest multiband superconductivity, with: i) a con- 
ventional phonon-mediated superconducting coupling (at 
least in non-magnetic ones); ii) weak intcrband coupling 
(in YNi2B2GiS'iS and in LuNi2B2Gi^ii^) with differ- 
ent gaps and T^, iii) a generally anisotropic distribution of 
gap values over the FS, with a larger gap along the c axis; 
iv) a separation of superconducting and magnetic order 
parameters on different sheets of the FS (in IIoNi2B2C^^^ 
and ErNi2B2Gi^). The research in this field is still going 
on and new measurements will certainly help clarifying 
this complex situation. 



E. PCAR in Graphite-intercalation compounds 

The recent discovery of the new superconduct- 
ing Graphite Intercalation Compounds (GICs) CaCg, 
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FIG. 28: (a) Low-temperature experimental PCAR spectra 
(circles) in CaCe for ab-plane (red) and c-axis (blue) current 
injection, superimposed to the relevant single-band BTK fit 
(dashed lines). Black solid lines are the theoretical conduc- 
tance curves at T = calculated from the full k distribution 
of gap values and Fermi velocities, and smeared with the indi- 
cated (experimental) F values, (b) Distribution of measured 
gap values and relevant gaussian fit. (c) Gaps vs. T extracted 
from the fit of the temperature dependence of the curves in 
(a). 



YbCe^, and SrCgi^ has renewed the interest for this 
long-known class of compounds. CaCg shows the high- 
est Tc = 11.5 K among this class of compounds. Its 
lattice is made up of alternating graphite layers and Ca 
planes, with a rhombohedral structure^i^. The similar- 
ity with MgB2 is striking and indeed the electronic band 
structurei^ includes the a and tt bands, though the for- 
mer are completely filled and play no role in supercon- 
ductivity. According to first-principle calculations i^iii^ 
superconductivity arises from the coupling of both Ca 
and C phonon modes to the carriers of the so-called in- 
terlayer band, which is formed by C and Ca orbitalsi^S. 
The Fermi surfacei2£ consists of 7r-band warped cylin- 
ders parallel to the c axis, and of interlayer-band FS 
sheets created by the intersection of some of these cylin- 
ders with the nearly-spherical Ca orbitals'^i'^ . Unlike in 
MgB2 (see FigfT^a)) the calculated gap^^ is continuously 
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distributed in energy between 1.1 and 2.3 meV and the 
mixing of C and Ca states prevents a true multigap su- 
perconductivity. However, the gap mapping on the FS 
shows that it changes from one FS sheet to the other and 
is, on average, sHghtly larger in the 3D interlayer band 
(actually depending also on the wavevector k within each 
sheet). This suggested us to use directional soft PCAR in 
CaCe single crystals to observe the predicted anisotropy. 
The extreme sensitivity of CaCe to air and moisture re- 
quired to cleave the sample and make the contact in inert 
atmosphere, and to seal the whole sample holder before 
transferring it to the cryostat. Due to the small Tc = 11.5 
K and the smallness of the effect to be observed, part of 
the measurements were carried out at 400 mK??. 

Fig lMl a) shows a c-axis and a a6-plane PCAR spec- 
trum at T = 0.4 K (symbols). The position of the con- 
ductance peaks is clearly different and indeed the fit with 
a single-band 2D BTK model gives Aaf, = 1.44 meV and 
Ac — 1.7 meV. The values of Z are systematically higher 
in c-axis contacts (0.74 < Z < 1.01) than in a5-plane ones 
(0.48 < Z < 0.75), in agreement with the different Fermi 
velocities in the two directions {vab = 0.54 • 10^ m/s, 
Vc ~ 0.29 • 10^ m/s). These spectra are in excellent agre- 
ment with ab-initio calculations of the Andreev conduc- 
tance carried out^ as described in Sect IV 71 eg 1311 To al- 
low a comparison with experiment, the theoretical curves 
at T = (shown in Fig. ^ were smeared with the experi- 
mental values of F (i.e. 0.6 meV for the a6-plane contact, 
0.8 meV for the c-axis one) neglecting the (much smaller) 
thermal smearing. The results are shown in Fig |28l as 
solid black lines. FigESTb) shows that the distributions 
of low-temperature gap values obtained in several c-axis 
and a6-plane contacts are approximately Gaussian and 
overlap only slightly. The temperature dependence of 
the gaps is shown, for the spectra in Figl^BTa). in panel 
(c). The gap always follow a BCS-like trend with gap 
ratio 2Aab/kBTc = 2.98 and 2Ac/kBTc = 3.48. All these 
results gave the first direct evidence of gap anisotropy 
in CaCe and showed a very nice example of a success- 
ful feedback between PCAR experiments and theoretical 
predictions from ab-initio calculations. 

VIII. CONCLUSIONS 

Multiband superconductivity was theoretically investi- 
gated since the Sixties, but it was often considered as an 
exotic although interesting possibility, with little practi- 
cal relevance, because of the feebleness of its effects even 
in the very few cases when they were detected, or the 
very small critical temperature of the materials that dis- 
played them'*^. On the other hand, point-contact spec- 
troscopy was invented in the mid-Seventies and, although 
initially applied to the study of normal metals, it was 
soon understood to be a powerful spectroscopic tool in 
the study of superconductors. However, other techniques 



such as STM and ARPES with better spatial and mo- 
mentum resolution largely dominated, in the Nineties, 
the study of cuprates. With the discovery of MgB2 in 
2001, multiband superconductivity suddenly became a 
promising and popular field of research, and PCAR spec- 
troscopy rapidly acquired a great relevance thanks to its 
quick and successful application to this compound. The 
amount of information this technique has been able to 
provide, even in samples that were impossible to analyze 
by STM and ARPES, was probably a surprise for most of 
the superconducting community. Since then, and thanks 
to its reliability, simplicity and flexibility, PCAR spec- 
troscopy has played an important role in the investiga- 
tion of the superconducting properties of many new (and 
less new) compounds. 

In this review we have tried to present the recent appli- 
cations of PCS and PCAR spectroscopy to the study of 
multiband superconductors, starting from MgB2 (either 
pure or doped) to continue with borocarbides, graphite- 
intercalation compounds and the recently discovered Fe- 
based superconductors. We have shown that PCAR mea- 
surements can provide information on the number, the 
amplitude and the symmetry of the superconducting or- 
der parameter(s), but also - when integrated with first- 
principle calculations and with the Eliashberg theory ~ on 
the coupling strengths (both within and between bands), 
the scattering rates, the densities of states, the k depen- 
dence of the gap within the various sheets of the Fermi 
surface, and so on. We have provided the reader with a 
simple theoretical introduction to PCS and PCAR spec- 
troscopy, showing how the limitations of the original, pi- 
oneering Blonder-Tinkham-Klapwijk model (often used 
even today to fit the point-contact spectra) can be over- 
come to improve the degree of approximation to the real 
case and to make PCAR a much sharper tool for the 
investigation of unconventional superconductors. 
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If the effective electron mass can be supposed to be the 
same for metals 1 and 2, the Fermi velocities in equation 
111! can be replaced by the Fermi wavevectors. 
In this calculation, interference effects between bands 
were not taken into account. In a recent papei«22., it has 
instead been shown that such effects can in principle give 
rise to observable features in the Andreev conductance 
spectra not only in iron pnictides, where the order pa- 
rameter changes sign on different bands, but also in MgB2 
where the order parameter has the same sign on both a 
and TT bands. 



